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| less than they could he.** J. Ditngcy fl 906P ' 

"Most of the important phenomena involve simultaneous 
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illllllllly; to make simultaneous measurements at ^veftd 
hundred locations." CE. A Mhvsin 11907P' 
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(n from “Inaugural Lecture as Professor of Physics at Imperial College,” 1966; quote courtesy of W. J. Plughes. 
<2) from “Comments and Speculations Concerning the Radiation Belts,” page 303 (Proceedings of the Joint 
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Executive Summary 

Magnetospheric Constellation Dynamic Response and Cou- 
pling Observatory (DRACO) is the Solar Terrestrial Probe 
(STP) designed to understand the nonlinear dynamics, re- 
sponses, and connections within the Earths structured 
magnetotail, using a constellation of - 50-100 distributed 
vector measurement spacecraft. DRACO will reveal 
magnetotail processes operating within a domain extend- 
ing 20 Earth radii (R E ) across the tail and 40 R E down the 
tail, on spatial and time scales accessible to global circula- 
tion models, i.e., - 2 R E and 10 seconds. Visualizing and 
understanding them will require observations analogous 
to a network of weather stations distributed within this 
domain. DRACO will reveal simultaneously for the first 
time both the global spatial structures and the time varia- 
tions of the magnetotail. It will determine which phenom- 
ena are responses to solar wind inputs and which occur as 
a result of internal instabilities. In particular, it will re- 
veal the locations and extents of the instabilities tha_t trig- 
ger the explosive release of solar wind energy, mass, and 
momentum stored within the magnetotail, how these en- 
tities are transported, and the means by which magneto- 
tail phenomena are propagated between regions and to the 
auroral ionosphere. 

The magnetotail is a critical volume of the geospace envi- 
ronment wherein global circulation of magnetic fields and 
plasmas is regulated in response to changing solar wind 
conditions. In it, impulsive localized flow bursts launch and 
dissipate, powerful electrical currents form and evolve 
abruptly, and magnetic energy is explosively converted to 
particle energy. The fundamental plasma process known 
as magnetic reconnection is thought to occur during 
substorms, an important building block of "space weather.” 
These are recurrent energy releases that become more fre- 
quent during magnetospheric storms. The dynamism and 
turbulent evolution of the magnetotail have humbled our 
efforts to observe and understand it using individual space- 
craft. The magnetotail magnetic fields and plasmas do pos- 
sess an underlying, slowly varying coherent structure, but 
strongly turbulent flow T s and fields are usually large com- 
pared to the mean field or flow. Thus, globally coherent 
pictures of the system dynamics become lost in the "noise” 
of individual measurements. Despite over 30 years of re- 
search with ever more sophisticated instrumentation on 
ever larger and more complex spacecraft, fundamental 
questions concerning the dynamic response of the 
magnetotail remain unanswerable. Accordingly, scientific 
progress has slowed. Intelligent, reasonable scientists can- 
not reach consensus on these issues, not for lack of models 
and theories, but because of a lack of relevant measure- 
ments. Neither current single spacecraft nor tight groups 
of spacecraft can resolve these fundamental controversies. 

To provide the first ever global time-evolving vector field 
and streamline images of this important region, DRACO 
will use rapidly developing technologies to deploy a "con- 
stellation” of nanospacecraft. With resources of ~ 10-20 kg 
and 10W apiece, 50-100 nanosatellites will be deployed in 
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highly elliptical, equatorial orbits with common perigees 
of 3 R f . and apogees distributed from 7 to 40 R E , yielding 
mean inter-spacecraft separation of - 1-2 R E . 

The nominal DRACO mission has a design lifetime of 2 
years, and is currently scheduled for a 2010 launch. While 
the enabling technologies for DRACO are developed or in 
development for the ST-5 Nanosat Constellation Trailblazer 
mission of the New Millennium Program, additional re- 
sources are required to demonstrate mass manufacturabil- 
ity and to optimize miniaturization while preserving func- 
tionality. With such resources, DRACO development could 
be accelerated with the launch as early as 2008. 

Magnetospheric Constellation DRACO represents the logi- 
cal outgrowth of a sequence of STP missions, designed to 
explore plasma transport and energy conversion processes 
over spatial sizes ranging from the distance to the Sun to 
the size of low-energy particle gyro-orbits. The Magneto- 
spheric Multiscale (MMS) mission will focus on the small- 
est scale, microphysical processes occurring within and 
near magnetospheric boundary layers: magnetic 
reconnection, charged particle acceleration, and eddy tur- 
bulence. It will serve as the plasma physical "microscope,” 
investigating scales too small to be resolved by global cir- 
culation models. In contrast, DRACO is designed to be a 
"meso/macro-scope” for the magnetotail. It will resolve per- 
sistent controversies by providing the required observa- 
tions. Ultimately it will yield a new understanding on 
which we shall build a predictive science of next-genera- 
tion magnetospheric meteorology, adding to our collective 
body of knowledge relating to complex nonlinear dynami- 
cal systems. 

Magnetospheric Constellation DRACO Highlights: 

* DRACO will answer the following questions: How does 
the magnetotail control energy flow? What processes 
control magnetotail structure and dynamics? How do 
the physical processes and regions couple over the hi- 
erarchy of scales? 

* Leads to closure of crucial magnetotail controversies 
that have remained unresolved for more than 30 years, 
owing to lack of coordinated multipoint observations. 

• Single Delta launch creates a magnetotail constella- 
tion with possible launch in mid-2010 or earlier. 

* 50-100 spacecraft with 10 sec time resolution, at mean 
separation of - 1-2 R £ over - 20 x - 30 R E domain, re- 
solve all features of global circulation models. 

• Primary science accomplished annually when the con- 
stellation sweeps through the magnetotail. Ancillary 
magnetospheri c/magnetosheath/solar wind science ac- 
complished during balance of each year. 

• New technology requirement is for miniaturization and 
mass manufacturability of nanosatellites and their in 
strument payloads. 
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1.0 Introduction and Science Objectives 

1.1 Magnetotail Configuration 

The geomagnetic field carves out a three-dimensional cav- 
ity in the flow of solar wind plasma, known as the mag- 
netosphere, as illustrated in Figure 1. Supersonic solar 
wind plasmas stream away from the Sun and are heated, 
deflected, and decelerated at the bow shock. Within the 
magnetosheath, the slowed plasma flows around the mag- 
netopause, the locus of positions where the incident solar 
wind pressure balances that of the geomagnetic field. The 
magnetopause is not a perfectly impenetrable boundary, 
though. A small fraction of the solar wind mass, energy, 
and momentum enter the cavity, resulting in the transfer 
of energy to the magnetosphere. The amount of energy and 
the locations for entry are controlled by constantly vary- 
ing solar wind conditions. The transferred energy is stored 
principally in the region called the magnetotail, avast re- 


gion of dynamically evolving magnetic fields and plasmas 
stretching far from the Earth in the anti-solar direction. 

The Earth's magnetotail consists of several distinct regions 
(ref. Figure 1). Much of the plasma in the magnetotail 
resides in the plasma sheet, a region of hot plasma ex- 
tending across the tail from the dawn to the dusk edge, 
and extending from the distant tail into the inner mag- 
netosphere. The thermal pressure of the plasma sheet is 
substantial and supports a delicate balance between the 
plasma sheet and the stretched field lines of the low-den- 
sity northern and southern tail lobes. Enshrouding the 
northern and southern lobes are the plasma mantles or 
High Latitude Boundary Layers (HLBL). The HLBL are 
regions of dense, cool plasma flowing along the magnetic 
field lines in the anti-Sunward direction. They contain 
plasma mainly of solar wind origin with a minor compo- 
nent that has escaped from the ionosphere. On the dawn 
and dusk flanks of the magnetosphere are the Low-Lati- 
tude Boundary Layers (LLBL), regions of anti-Sunward 



Figure 1 . The figure illustrates the computed magnetotail and Its major component regions and boundaries in three 
dimensions; figure was generated from the results of magnetohydrodynamic global simulation modeling. Color-coding 
indicates plasma density on intersecting planes at Y=0, and at X '= 40 Earth radii. Labels indicate significant bound- 
aries and other features. Light blue magnetic field lines originate from geomagnetic mid -latitudes, while dark blue 
field lines originate in the polar cap region. Red dots indicate typical constellation spacecraft positions near the plasma 
sheet (green color on the Y=0 plane). Figure is provided courtesy of J. Raeder, UCLA. 
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flowing plasma of solar wind origin, again mixed with mag- 
netospheric plasma. Finally, the northern and southern 
Plasma-Sheet Boundary Layers (PSBL) lie between the 
plasma sheet and the lobes. The PSBL typically contains 
magnetic-field-aligned beams of electrons and ions that 
originate from magnetic reconnection farther down the tail. 
These seemingly disparate regions are globally intercon- 
nected. This intricately linked system responds dynami- 
cally to external boundary conditions as well as to inter- 
nal reconfigurations on short time scales. 

Much of the plasma of mixed internal and external origin 
becomes trapped within the geomagnetic field and is un- 
able to escape directly downstream. This plasma partici- 
pates in the global circulation of the magnetosphere. In 
the closed field region of the magnetotail, the magnetic 
field lines are stretched into a long plasma tail. Curvature 
forces in the stretched field lines ultimately cause the en- 
trained plasmas to convect Sunward back toward the day- 
side magnetosphere from which they came, thereby com- 
pleting the cycle of global circulation. The earliest models 
of the magnetosphere cast this process as a simple time- 
independent circulation. Today, we know that the Earth- 
ward, return flow is both highly variable in time and highly 
structured in space. Consequently single spacecraft stud- 
ies have faced overwhelming challenges in unambiguously 
describing the fundamentals of this flow. The flow is a vi- 
tal link in understanding how the solar wind energy propa- 
gates through and thereby affects geospace. Accordingly 
it is a focal point of the DRACO mission. 

With the exception of global magnetohydrodynamic (MHD) 
models, global magnetospheric models are static 
( Tsyganenko, 1989, 1995]. The large magnetic field vari- 
ability in the databases that are used as input to those 
models [Stern and Tsyganenko, 1992] suggests that a sig- 
nificant part of the physics that determines the most prob- 
able magnetospheric state is missed. Most often, severe 
modification of the model currents is necessary away from 
expected state [i.e., for the measured Kp or AE parameter] 
when actual time-dependent situations are considered 
[ Pulkkinen , 1991; Baker et al,. 1993]. It is impossible for a 
static model parameterized on global Indices to capture 


the complexity the richness and the physics of the instan- 
taneous magnetospheric configuration. However, a large 
array of measurements from magnetospheric probes can 
produce, when inverted using techniques available from 
statistical [Tsyganenko and Usmanov, 1982], simulation 
\Berchem et aL, 1995] and atmospheric [Ghil and 
Ma lanotte -Rlzzoll, 1997] models, an instantaneous "image 7 ' 
of currents and fields that is consistent with the average 
tail structure but representative of the instantaneous mea- 
sured fields. This can produce a far clearer picture of the 
magnetospheric current systems, their sources and sinks 
and their temporal evolution. Parameterization versus 
substorm-phases and upstream conditions can produce the 
"typical” substorm sequence, i.e., a time history of the mag- 
netospheric system's response to qualitatively different 
drivers. 

1.2 The ‘Quiet’ Plasma Sheet 

Global Versus Turbulent Transport and Thermody 
namics — Magnetospheric Constellation DRACO will pro- 
vide the critical new observations for visualizing and un- 
derstanding large-scale plasma transport within the mag- 
netotail. Single spacecraft missions have explored this re- 
gion comprehensively over time, but we still have only a 
poor understanding of its global, quiet-time, instantaneous 
behavior. Indeed, past single-satellite observations have 
presented important glimpses into the underlying texture 
of the near-Earth magnetotail. It is known, for instance, 
that the magnetic field of the plasma sheet has large fluc- 
tuations and that its plasma flows are highly erratic 
[Hayakawa et al, 1982; Borovsky et al, 1997]. 

While the most dramatic flow events are associated w'ith 
geomagnetically active times, the most probable state of 
the plasma sheet (90-95% of the time) is the state possess- 
ing no significant bulk flow [Baumjohann et al, 1989, 1990]. 
That state exhibits random flow directions with peak-to- 
peak amplitude many times larger than the average flow 
itself [Angelopoulos etal, 1 993]. Despite the large flow vari- 
ability, the average flow pattern exhibits behavior that can 
be understood in simple 2-fiuid theory [Angelopoulos et al., 
1993] or single particle orbits [Spence et al, 1993a, 1993b] 
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IRM DATA (Angelopoulos, 1996) 
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Figure 2. The mean flow and variability of the magnetotail showing that a Sunward /duskward circulation pre- 
vails (in the positive X and Y directions), but that the variability of the flow is as large or larger than the mean flow 
speeds (right panel). Figure is courtesy ofV. Angelopoulos. 
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(Figure 2). Significant variability about the magnetic field 
average also exists [Stern and Tsygancnko, 1992], render- 
ing orbit integration schemes vulnerable to criticism 
\CatteII et al, 1995]. Is the variability the consequence of 
nearby high-speed flows [ Angelopoulos et al. , 1994; 
Borovsky et a I., 1997] or the result of nonadiabaticion tra- 
jectories under steady external conditions [Ashour-Abdalla 
et al, 1995]? Both theories are viable yet presently 
unresolvable with available observations. DRACO will pro- 
vide the comprehensive view of plasma sheet flows needed 
to unravel these questions. Furthermore, what role does 
the flow variability have on resolving the long-standing, 
"pressure balance inconsistency" which arises in models 
of time-independent, laminar tail convection? 

As noted above, statistical treatments reveal an average 
Earthward flow in the plasma sheet, but with fluctuations 
as large or larger than the average flow. Indeed, this ir- 
regular flow of the plasma sheet may become fully turbu- 
lent. Turbulence data-analysis methods have been applied 
to single-satellite measurements of fields and flows 
[Borovsky et al, 1997]. Models of the magnetotail morphol- 
ogy based upon turbulent transport have been developed, 
and models of auroral activity driven by flow turbulence 
in the plasma sheet have been constructed [Swift, 1981; 
Song and Lysak, 1988]. 

Velocity auto-correlation times are of the order of 2 to 3 
minutes and mixing lengths of the order of 1-2 R £ . This 
suggests that although large-scale eddies are occasionally 
observed in the plasma sheet [ Hones et al, 1978, 1981] 
underlying effective mixing also takes place at shorter scale 
lengths. In order to characterize this state of the plasma 
sheet, correlative measurements at many points of the sys- 
tem are necessary in order to measure the wavelength spec- 
trum of the plasma sheet flows. This will tell us whether 
substorm processes generate the variability or whether it 
is an inherent, fundamental property of the convecting 
plasma sheet. Distribution functions are necessary in or- 
der to ascertain whether the flow variability is due to par- 
ticle bunching or a velocity shift of the entire particle 
distribution. 

Recent results from the International Solar Terrestrial 
Project (ISTP)-era indicate that beyond - 30 R £ the mag- 
netic field configuration is tail-like and the plasma beta is 
high. Present observations indicate that in this region the 
flow is reasonably laminar. Closer in, the Earthward-jet- 
ting plasma encounters the deflecting effect of the Earth s 
rigid dipolar field [ Baumjohann , 2000]. In this transition 
region, the flow is observed to grow more turbulent. The 
transition from laminar to turbulent, which may addition- 
ally involve bursty and localized reconnection, can only be 
realistically studied using a constellation of spacecraft. 

The fundamental thermodynamic properties of the plasma 
sheet are also not clearly understood. Huang et al. [1992] 
have shown that plasma sheet temperature increases af- 
ter substorm onset represent an important plasma sheet 
response to geomagnetic activity but that they do not al- 


ways correlate with flow enhancements. Although single 
particle approaches [Spence, 1993a, 1993b; Ashour-Abdalla 
et al, 1994] have led to certain predictive pictures of the 
plasma sheet pressure, temperature and density spatial 
profiles, such profiles have not been reported except in en- 
ergetic particles [Krimigis and Sarris, 1979]. Thus, there 
is a need to observe the system and its evolution at many 
points simultaneously for a given thermodynamic state. 
This entails fine spatial resolution observations of the in- 
stantaneous flow pattern, density and temperature to ob- 
serve the source and propagation of the resulting heating. 

The reason for the fine spatial resolution is twofold: First, 
the acceleration regions are localized in space [Krimigis 
and Sarris, 1979; Sergeev et al, 1986b, 1996; Angelopoulos 
et al, 1997]. Second, the plasma sheet flow and thus par- 
ticle trajectories are variable at inferred spatial scales of 
- 2 R e . Without a dedicated constellation mission, it has 
been shown to be unlikely that two spacecraft will have 
many chances of crossing the same flow line at a given 
time. The high spatial density and high-time resolution 
distribution functions provided by DRACO are necessary 
to resolve this problem. 

Attempts to characterize the detailed thermodynamics of 
the plasma sheet (such as measuring its polytropic index) 
have been made [Baumjohann and Paschrnann, 1989; 
Huang et al, 1989; Coertz and Baumjohann, 1991]. How- 
ever, there is always a great uncertainty in any inferred 
result, because with a single spacecraft we cannot obseive 
the evolution of a single flux tube. Multiprobe studies that 
utilize two (or more) spacecraft traversing plasma stream- 
lines will be able to answer the aforementioned question. 
However, discoveries from such a multiprobe mission may 
result in a complete change of mindset regarding this is- 
sue; if the plasma flow is nonlaminar, the concept of a 
streamline may be invalid. 

1.3 Global Aspects of Magnetic 
Reconnection 

Magnetic reconnection is one of the most important en- 
ergy transfer mechanisms operating throughout the mag- 
netosphere, yet we have little knowledge of how it is trig- 
gered and what its consequences are globally. It is thought 
to be the driver by which energy enters the magnetosphere 
from the solar wind, and important for the ultimate re- 
lease of this transferred energy which builds up episodi- 
cally in the magnetotail. Over long enough time intervals 
and integrated over the entire spatial volume, dayside and 
nightside reconnection rates must balance so that no mag- 
netic flux is permanently added to or removed from the 
Earths magnetosphere. Transient increases of one rate 
over the other are important drivers of short-term magne- 
totail dynamics in which the magnetotail moves from one 
state to another. The magnetotail varies in response to con- 
stantly varying solar wind conditions imposed on the mag- 
netopause boundary. Some of these responses occur in di- 
rect relation to the external drivers, while others are tied 
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to internal global instabilities. Despite the vastness of the 
three-dimensional magnetotail, its regions are closely 
linked electrodynamically, and respond globally in complex 
ways. The Magnetospheric Constellation DRACO mission 
is designed to resolve these important features of the 
magnetotail. 

Based on the analysis of a vast quantity of single space- 
craft data, the following is a phenomenological model of 
how the solar wind driver, through magnetic reconnection, 
sets the magnetotail toward a global instability known as 
a magnetospheric substorm. Magnetic reconnection at the 
dayside magnetopause breaks open closed geomagnetic 
field lines and "reconnects” each end to solar wind field 
lines, coupling and mixing the solar wind and internal plas- 
mas. The motion of the solar wind drags the plasma on 
these "open” field lines downstream and into the 
magnetotail. The field lines are added to the northern and 
southern lobes and stored there as excess magnetic en- 
ergy, further squeezing and stretching the fields in the 
plasma sheet. 

The stretching continues until some yet undetermined in- 
stability occurs. The excess open field lines in the lobe must 
eventually reconnect a second time to conserve flux; this 
time going from an open to a closed topology. It is this re- 
connection in the tail between magnetic flux tubes once 
concentrated in the lobes that most likely powers geomag- 


netic substorms — the most common dynamical 
magnetotail event. The solar wind and ionospheric ends of 
lobe field lines decouple, releasing stored energy, and re- 
storing closed field lines within the magnetosphere and 
interplanetary magnetic fields in the solar wind. The ex- 
tended closed magnetic field lines, and the plasma on them, 
then collapse back closer to a dipolar configuration. The 
process creates a magnetic bubble of plasma, a plasmoid, 
downstream of the reconnection region that moves back 
into the solar wind \f loshinoet a1. t 1996; Mukai et ah, 1996; 
Slavin et ah, 1998]. Figure 3 is a cartoon depicting the 
gross morphology of the magnetotail just after the global 
substorm instability has occurred. 

While the stretching, growth phase of a substorm can be 
gradual (~ 1 hour), the relaxation toward a dipolar con- 
figuration, or dipolarization, is typically abrupt and dra- 
matic (~ 1 minute), characterized by high inductive elec- 
tric fields and fast plasma flows. The magnetic energy once 
stored in the lobe fields is converted rapidly at substorm 
onset to particle energy that is then released through the 
plasma sheet. Some particles are accelerated to high ener- 
gies and injected into the inner magnetosphere. Others pre- 
cipitate to form auroras at high latitudes. 

While initially a regionalized phenomenon, magnetospheric 
substorm effects become demonstrably global and evolve 
rapidly through the entire volume. Syntheses of single point 
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Figure 3. A cartoon of the global consequences of tail reconnection during a substorm. After a 
reconnection site has formed near Earth, fast flows propagate Sunward accompanied by a rapid 
dipolarization. A plasmoid has formed and is traveling downstream , creating a traveling com 
pression region in the magnetic lobes. DRACO nanosatellites are shown schematically through 
out the plasma sheet. Courtesy of J. Slavin/NASA Goddard Space Flight Center (GSFC). 
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measurements have led to statistical cartoons such as that 
shown above. While significant insights into tail dynamics 
have come from such analyses, our views have been lim- 
ited and competing theories have remained untested. To 
compensate for this observational deficit, further insights 
into substorm dynamics have been gleaned by appealing 
to MHD models. The macroscopic structures illustrated 
schematically in Figure 3 are present in virtually all MHD 
models of the dynamic magnetotail (Figure 4). These nu- 
merical models reveal intricacies of the magnetic topolo- 
gies and complexities in the flows not only during substorm 
intervals, but also during the relatively quiet resting state 
of the magnetotail. They also predict dramatic variability 
of their time evolution that certainly confounds analysis 
by single, or even several, tightly clustered spacecraft. From 
the MHD model results, it is clear that magnetic recon- 
nection plays a major role in the configuration and dynam- 
ics of the magnetotail, and that reconnection is likely both 
bursty and patchy. 

The magnetospheric research community relies on a con- 
cept of magnetotail structure that has been constructed 
by averaging measurements taken by various single sat- 
ellites. In large measure, w 7 e rely on understanding the 
magnetotail's global dynamics through computer simula- 
tions having relatively coarse spatial resolution and some- 
times incomplete parameterization. The community has 
come to realize that deficiencies of this magnetotail 
conceptualization are impeding progress in understand- 
ing how the magnetosphere actually behaves. DRACO will, 


for the first time, reveal in three dimensions the true in- 
stantaneous evolution of the elusive substorm sequence in 
the magnetotail that we have only seen glimpses of in 
sparse data and physical models. 

1.4 The Substorm Problem: Chronology 
and Associated Mesoscale Structures 

As detailed earlier, the Earth’s magnetotail is a complex 
and dynamic region of the magnetosphere. Intense, small- 
scale plasma flows typically conceal the large-scale flow 
pattern withi n the plasma sheet. The magnetic field in the 
plasma sheet often exhibits small-scale kinks and tangles. 
However, during certain intervals, the Earth’s magneto- 
sphere undergoes large coherent changes in its global mor- 
phology, events we call substorms and magnetic storms. 
These phenomena are not fundamental physical processes 
in themselves, but rather a collection of chronologically or- 
dered physical processes that repeat episodically in some 
quasi-regular way. MHD models suggest that there is sub- 
stantial coherence within these global reconfigurations, a 
feature that is extremely hard to arrive at from single point 
measurements. The meta-process of a substorm represents 
the most frequent mode of significant magnetospheric dy- 
namics, and thus provides a practical context for under- 
standing the consequences of the underlying microphys- 
ics. While substorms have galvanized the community’s in- 
terest, a consensus understanding of the substorm has been 
hampered by data sparseness. 
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Figure 4. A snapshot from an MHD simulation of the dynamic magnetosphere. 
Sunward of a reconnection site near Earth, field lines are rapidly depolarizing , while 
a recently formed plasmoid is being ejected farther down the tail. Courtesy of J Raeder/ 
UCLA. 
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1,4.1 Chronology and Causality 

Ever since Syun Akasofu [1964] showed that auroral phe- 
nomena could be organized into characteristic recurrent 
patterns called "auroral substorms,” substorms have been 
a topic of vigorous debate. In the ionosphere, a global pic- 
ture can be obtained with all-sky cameras, radars, and mag- 
netometer networks; in space, we depend on statistical 
analyses to compose a time history of events. There is broad 
consensus that the magnetotail exhibits a characteristic 
behavior during magnetospheric substorms, defined by 
diverse observed aspects including auroral expansion, mag- 
netic field dipolarization, current wedge formation, 
dispersionless particle injections, bursty bulk flows, and 
current disruption. There is, however, no consensus on the 
physics of substorms, as suggested by the persistence of a 
designation that is only useful for classification and con- 
veys little understanding of physical basis. 

Recently, the debate has centered on the first few minutes 
of the substorm expansion phase: the so-called "substorm 
onset.” The first sign of substorm activity should be located 
close to the causative source. From a few fortuitous 
multispacccraft conjunctions, we have learned that the field 
dipolarization starts at a localized region of the order of 1 - 
2 R e \Ohtani et ah, 1992]. The field dipolarization propa- 
gates longitudinally \Nagai, 1982] at a speed of an hour of 
local time every 1-3 minutes and radially tailward \Jacqucy 
etaL, 1993; Ohtanietal , 1991, 1992] with a speed of- 250 
km/s, or 2.5 RJminute. Tailward flows are also observed 
during substorms and most often at distances greater than 
19 R,.; Earthward flows are observed closer to Earth 
\Baumjohann etaL , 1990; Angelopoulos etaL, 1994]. These 
flows have been interpreted as signatures of magnetic 
reconnection. However, the time and location of the first 
magnetospheric signature of substorm onset and its rela- 
tionship with the ground observables still remain a mys- 
tery. The present lack of an observational resolution of this 
point, in turn, inhibits theoretical progress in the under- 
standing of the physics that govern the phenomenon. Two 
theories, each associated with a different initial location of 
substorm onset have been proposed, and each predicts a 
different chronology and causality that is testable by the 
DRACO mission. 

One theory [Lui, 1991, 1996] holds that substorms start 
near or slightly outside geosynchronous orbit. Theories in 
this class suggest that the electrical current that produces 
the stretching of tail field lines is disrupted from flowing 
across the tail, instead flowing down along magnetic field 
lines and through the ionosphere. This forms the "substorm 
current wedge” [ McPherron et a]., 1973] — a transient 
substorm feature inferred from magnetometer measure- 
ments near geosynchronous orbit. The new current sys- 
tem allows the local magnetic field to relax and become 
more dipolar, launching a rarefaction wave that propagates 
down the magnetotail [Jacqiwy et al. , 1993; Lopez et aL, 
1993]. Some 20-30 Earth radii down the tail, it initiates 
increased magnetic reconnection of the type that turns open 


lobe field lines back into closed field lines in the magneto- 
sphere. This causes plasma to rush Earthward through 
the plasma sheet, allowing the plasma sheet to quickly 
reestablish itself and recover its plasma pressure and cross- 
tail current. 

The other theory [e.g., Schindler, 1974; Birn et al., 1996; 
Baker et al., 1996; McPherron et al., 1973; Hones and 
Schindler, 1979] argues that substorms start with increases 
in reconnection at a distance of 20-30 R t . The stretched 
lobe magnetic field is reconnected, creating a burst of 
plasma flow that launches captured plasma and field to- 
ward the Earth, at the same time forming a plasmoid that 
is shed downstream \Slavin etaL, 1993; 1998]. The Earth- 
ward burst launches a compression wave toward Earth 
[Moore et al., 1981 \ Reeves et al., 1996 ,LietaL, 1998], cre- 
ating a flow of particles with a wide range of energies (100 s 
of eV to 100s of keV) that runs into denser plasma and 
stronger magnetic field and is slowed and deflected as a 
result. The braking of the flow creates the observed 
substorm current wedge in this scenario. Also, at the same 
time, energetic particles up to - 500 keV are injected from 
the inner edge of the magnetotail to the trapping regions 
of the inner magnetosphere. 

Both of these theories seek to account for the characteris- 
tics of a substorm in the magnetotail, including the 
dipolarization (relaxation) of the magnetic field to a less 
stressed state, the substorm current wedge, and the accel- 
eration and injection of plasma into the dipolar region. 
However, the two have completely opposite flows of cau- 
sality: one case, from the tail toward the Earth; the other, 
from near the Earth out into the magnetotail. Figure 5 
shows these two models schematically [Baumjohann, 2000] . 
The difficulty lies in the localization of the initial phenom- 
ena both in space (1-2 R E ) and in time (10-20 sec) and its 
subsequent propagation over great distances (10-20 R g ). 
Single spacecraft within the vast magnetospheric system 
have little chance of being at the key region at the exact 
time of substorm onset. More importantly, the lack of si- 
multaneous measurements from adjacent "pixels” of this 
vast system prohibits placing the observed phenomenon 
in the context of the global magnetotail evolution. Many 
simulation results support the idea that the initial 
substorm disturbance originates farther out in the mag- 
netotail and then propagates inward as described above 
in the second theory. Limited observations [Ohtani, 1998] 
show that the field configuration changes propagate in- 
ward at a speed less than the sound speed, suggesting that 
they are not accomplished by a hydromagnetic compres- 
sion wave. Other observations suggest that the changes 
are simply convected by fast plasma flows (100s km/s), 
but field configuration accompanied by energetic particle 
injections were observed to travel as slowly as 25 km/s from 
6.6 R e to 4.5 R E . Lack of consistent observational confir- 
mation indicates that neither theory can be declared false, 
nor accepted as definitive. Thus, the most prevalent and 
important physical process in the magnetotail continues 
to escape scientific identification. 
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Figure 5. DRACO will resolve a long standing ambiguity in substorm chronology, illustrated in this figure. It will determine whether 
the magnetotail reconfigures rarefactively, from the inside out (right panel) or compress ively, from the outside in (left panel). 


The maximum altitude of the spacecraft necessary to cap- 
ture the reconnection process is determined from previous 
statistical studies of fast plasma sheet flows. Such studies 
\Hayakawa eta I., 1982; Nakamura etal, 1994; Hones and 
Schindler, 1979] indicate that at a distance of - 35 R E any- 
where between 25% and 80% of fast flows observed were 
tailward; the variability is due to the selection criteria and 
methodology. Recent Geotail results place the most likely 
position of observation of the near-Earth neutral line at 
around X=-25 R £ [Nishida etal., 1995] while they confirm 
that on the average the flows are tailward at distances of 
30-60 R e [Frank and Paterson, 1994]. Since tailward mo- 
tion of the neutral line occurs, the neutral line is likely to 
start most often Earthward of X=-40 R £ at substorm on- 
set. Thus, it is in the region between 7-40 R £ that onset and 
evolution of current disruption and near-Earth neutral line 
can be compared during the time of expansion phase onset 
of magnetospheric substorms. A suite of spacecraft at cross- 
tail separations ranging between 1 and 20 R & and at downfall 
position between 7 to 40 R £ are a necessary and sufficient 
condition for the study of the magnetospheric evolution dur- 
ing the course of a substorm. Quantities necessary and suf- 
ficient are high time resolution (1-5 s) magnetic field and 3- 
D ion moment measurements. 

1*4.2 Mesoscale Processes and Structures 

Particle Acceleration Fronts — Substorms are a primary 
mechanism of energization of the nightside ionosphere and 
the ring current. Although the ultimate source of energy 
during substorms is the solar wind, the transfer of that 
energy to the plasma sheet particles is not well understood. 
Modeling of geosynchronous injections during substorms 
as broad spatial fronts of Earthward-collapsing particles 
[Moore et ah, 1981; Mauk , 1986] (based on measured elec- 
tric fields observed in the near-Earth plasma sheet [Aggson 
etal, 1983; Maynard etal ., 1996]) has produced fairly good 
agreement between theory [Quinn and Southwood, 1982] 
and observations. In that scenario, it is the braking of the 
Earthward flows \Haerendel, 1992] resulting from either 
the current disruption or the reconnection process that 
causes the near-Earth signatures of substorms. Until now, 
flow braking is established only statistically in that the 
occurrence rate drops with distance along the tail. DRACO 
can establish the instantaneous pattern of this phenom- 
enon, which is a vital clue to its source. In the flow-brak- 


ing region (10-15 R £ ) P the flow causes pressure gradients 
and, perhaps, vorticity, which is the source of the Substorm 
Current Wedge (SCW). The azimuthal gradient has been 
established statistically by observations and in simulations. 
DRACO will for the first time establish the instantaneous 
pattern. 

Current Wedge — The longitudinal spread of the substorm 
current wedge [ Nagai , 1982] is coupled to the evolution of 
the westward traveling surge in the ionosphere [Roux et 
al., 1991; Robert et a!., 1984; Maynard et ah, 1996]. The 
extent in Y of the SCW at distances far from geosynchro- 
nous is neither known, nor possible to document with single 
spacecraft because it is a 3-dimensional system. Lopez and 
If// [1990] used AMPTE/IRM and Charge Composition Ex- 
plorer (CCE) observations to show that the SCW is irregu- 
lar in space. In addition, substorms take place most often 
as a series of consecutive, isolated activations. The near- 
Earth and ground response is the integral response to the 
resulting consecutive elemental current wedges. Such in- 
dividual activations have been proposed to be the elemen- 
tal processes [Lui etal., 1991; Sergeev et al . , 1986a, 1986b] 
that compose substorms. Multispacecraft observations will 
result in resolving the ambiguities associated with the 
number, location and extent of such activations. 

In fact, ISTP studies [Angelopoulos et ah, 1997; Slavin et 
al., 1997; Fairfield et al., 1998] are already providing di- 
rect measurements of the extreme localization of the el- 
emental activations and their relationship to the iono- 
sphere. Despite the limited number of near-by conjunctions 
of ISTP spacecraft while in the magnetotail, it is possible 
with these to glimpse the complexity and richness of the 
magnetotail plasma when studied at intermediate spatial 
scales. Perhaps, though, the most important lesson learned 
from ISTP is the scarcity of fortuitous conjunctions and 
the need for a distributed satellite system even when the 
goal is the study of the microstructures and macrostruc.- 
tures. The subliminal, and prevailing, understanding is that 
the current wedge maps to the entire width of the aurora 
and that its outermost boundary is the location where the 
current feeds from the magnetosphere to the ionosphere. 
Thus, the SCW is an important key in the chain of energy 
from the plasma sheet plasma to the ionospheric Joule 
heating. Understanding the three-dimensional evolution 
of the substorm current wedge system promises to solidify 
estimates of the magnetospheric-ionospheric energy cou- 


8 


Magnetospheric Constellation Mission... 


pling during a substorm, and in addition probe simulta- 
neously its component activations and their ionospheric 
mapping. 

“Bursty Bulk Flows” — The most prominent plasma 
sheet transport phenomena occur during short intervals 
and correlate well with substorm activity. Short-lived, high- 
velocity “bursty bulk flows” (BBF's) are observed within 
| Y | <10 R e from the mid tail axis [Baumjohann et al, 1990; 
Angelopoulos et al, 1994]. In the magnetotail tailward of 
X=-15 R e , and close to the neutral sheet, fast plasma sheet 
flows are thought to be a significant if not dominant part 
of the long-term average transport of the plasma sheet 
energy, particle and magnetic flux \ Angelopoulos et ah, 
1 994]. Present estimates for the important spatial and tem- 
poral scales are about - 2 R E and - 1 sec for flow bursts. 

A significant impediment to understanding their signifi- 
cance is the uncertainty' in their physical dimensions. With- 
out knowledge of that critical parameter, we have only a 
vague understanding of how much flux and plasma they 
transport. Are the fast flows that are observed at substorm 
onset of sufficient energy to account for substorm 
energization of the near-Earth environment and the iono- 
sphere? The answer depends on the 3-D evolution of the 
magnetotail flows. Given the extreme localization in Y of 
some of the fast flow events observed so far \Sergeev et al, 
1996; Angelopoulos et al, 1996; Angelopoulos etal, 1997], 
simultaneous multipoint measurements over an extent of 
dY -20 R e and at a probe density of 1 measurement per 1- 
3 R E in Y are necessary. 

MHD models demonstrate similar mesoscale features and 
an example is shown in Figure 6. Four panels are shown 
of the equatorial plasma sheet. Flow vectors (arrows) and 
parallel current densities (color coded) are shown through- 
out the magnetotail at four different times during a 
substorm simulation. A flow burst is seen to race Earth- 
ward in the post midnight region and then brake as it 
reaches the inner magnetosphere. As the flow brakes, a 
current wedge appears that couples to the ionosphere. 
While the effects of this flow burst and wedge are promi- 
nent in the post-midnight region, little effect is seen in the 
midnight or premidnight sectors. During this same run, 
other flow bursts (not shown) are seen to originate in the 
tail and move Earthward but always at a different loca- 
tion. Without the coverage of a constellation, the scope and 
significance of these flow channels remain unknown. 

Because the fast flows are bursty, they may likely come 
from localized reconnection. DRACO can determine the 
extent of the presumably local X-lines and will also be able 
to observe whether two or more small-scale X-lines can co- 
exist at any one time and whether these X-lines can de- 
velop into more global ones during major substorms. 
DRACO will provide the spacecraft spacing and distribu- 
tion necessary to unravel both the spatial extent and tem- 
poral evolution of BBF's as well as their detailed role to 
the substorm. 


Current Sheet Thickness and Evolution — Present 
theories that attempt to address the substorm onset ques- 
tion on the basis of current sheet instabilities require di- 
rect measurements of the current sheet thickness and in- 
tensity [e.g., Lui et al, 1991, Wang et al., 1990]. Attempts 
to model the current sheet thickness and density by in- 
verting the observed magnetic field measured at more than 
one location during fortuitous conjunctions of spacecraft 
has lead to some success [Pulkkinen et al, 1991]. In addi- 
tion, current measurement techniques in the tail using two 
spacecraft measurements [ McComas et al. , 1986] have 
given incontrovertible evidence for the existence of thin 
current sheets and their importance during substorms. 
Flowever, the only possible way to ensure that the tail-cur- 
rent system is adequately monitored at the time and place 
where it actually becomes unstable is to simultaneously 
monitor its thickness at many different downtail distances. 
This necessitates a CLUSTER-like system at many differ- 
ent distances. One way to achieve that is to design for some 
orbits with enough separation in Z (~ 1-2 RJ at different 
downtail distances. For current sheet densities of the or- 
der of 1 nA/m 2 which are typical of the cross-tail current, 
and intersatellite separations of 1 R E such measurements 
can be made comfortably with 0.5 nT accuracy, i.e., well 
within the capabilities of conventional magnetometer 
designs. 

Plasmoids and Flux Ropes — Plasmoids are large three- 
dimensional segments of the plasma sheet which are 
ejected down the tail in association with reconnection at 
the near-Earth reconnection site. When spacecraft are po- 
sitioned in the tail beyond X - -25 R F during magneto- 
spheric substorms, plasmoids are detected with great regu 
larity [e.g., Taguchi et al, 1998; Moldwin and Hughes, 1992; 
Slavin etal , 1993]. These large structures possess typical 
volumes and downtail speeds of - 1 0 3 — 1 0 4 cubic R E and 
- 400-800 km/s, respectively. These dimensions have been 
inferred from a synthesis of single spacecraft observations, 
but we lack a clear view of the actual origin, structure, and 
evolution of plasmoids. 

Schindler \ 1974] and Hones [1976] predicted the plasmoid 
phenomenon on theoretical grounds. As depicted earlier 
in Figure 2, these complex structures begin to form late 
in the substorm growth phase when reconnection at the 
near-Earth neutral line (NENL) is still limited to the closed 
magnetic flux tubes threading the plasma sheet [Taguchi 
et al ., 1998]. The location of the NENL near the time of 
substorm onset is now known to be X - -20 to -30 R E \Nagai 
et al, 1 998]. Plasmoid ejection begins at the moment when 
the last closed field line in the plasma sheet merges, thus 
allowing the lobes to come together and begin reconnect- 
ing "open’' field lines. It is the reconnection of the lobe field 
lines which liberates the large amounts of magnetic en- 
ergy necessary to power substorm expansion phase phe- 
nomena in the inner magnetosphere and the ejection of 
plasmoids down the tail. 
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Figure 6. Four snapshots from an MHD simulation of the dynamic magnetotail; frames are approximately 1 minute apart. A flow 
burst moves from the deep tail toward the Earth where the flow brakes. As the flow slows, a current wedge forms that feeds current 
to the ionosphere. Such regionalized features recur sporadically in the modeled magnetotail during substorms. Figure Is courtesy of 
C. Goodrich, UMd. 

The classic plasmoid model invoked a 2-dimensional ge- plasma data sets. The high-resolution measurements re- 
ometry, where the structure was simply a closed bubble of turned by Geotail clearly show that plasmoids grow rap- 


magnetic flux. The actual internal magnetic structure of 
plasmoids has been found to be quite variable; many pos- 
sess a helical structure and intense core magnetic fields 
resembling magnetic flux ropes [Moldwin and Hughes, 
1992; Slav in et al, 1998]. Such a spiraling of the plasmoid 
magnetic field had been predicted earlier from topological 
arguments \Hughes and Sibeck, 1987] involving the effects 
of the modest B y magnetic field component in the tail on 
reconnection. However, few details concerning the 3-dimen- 
sional structure of plasmoids and their magnetic connec- 
tivity to the rest of the tail or how they evolve with time 
are known. However, the overall impression was that 
plasmoids exhibited "compression regions" in the lobes 
\Slavin et al, 1993]. 

The simple-minded static view of plasmoid formation and 
ejection has recently been overturned by new ISTP-era re- 
sults. leda et al 1 1 998] has presented the first comprehen- 
sive survey of plasmoids in the Geotail magnetic field and 
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idly in all three dimensions following their release. A par- 
ticularly startling result is that Geotail frequently observes 
fast plasma flows within "young" plasmoids (i.e. t at X>- 
60R e ) with Vy - 100 to 200 km/s which produces a rapid 
plasmoid expansion in the east-west direction following 
its release. 

The Magnetospheric Multiscale mission with spacecraft 
separations in the tail of ~ 10 3 km will resolve the small- 
scale structure of plasmoid high field core regions. How- 
ever, measurements from a constellation of spacecraft at 
distances of X - 10 to 40 R E are necessary to determine the 
origins, dimensions, magnetic structure, stress balance and 
early evolution of plasmoids as well as their interaction 
with the surrounding magnetotail. The magnetic field, 
plasma and energetic particle measurements made across 
this region with spacings of 1-2 R £ will allow for determi- 
nation of the nature of the neutral lines, magnetic field 
reconfiguration and plasma sheet flow field which lead to 
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plasmoid formation and ejection. This constellation of 
spacecraft will also provide direct measurements of the 
plasma flow and circulation within the plasmoid, larger 
scale magnetic field topology, the rate of plasmoid expan- 
sion and any exchange of plasma with the external tail 
and, possibly, the magnetosheath. The combined magnetic 
and plasma data sets will enable direct examination of lo- 
cal stress balance between the regions within the plasmoid 
and between the plasmoid and the rest of the tail. Ener- 
getic ion observations will be used to map the separatrix 
layers associated with the NENL and Distant Neutral Line 
(DNL) and determine their location relative to the indi- 
vidual plasmoids. The energetic ion data will also be used 
to remotely sense the motion of boundaries/regions using 
well-known finite Larmor radius techniques [e.g., Owen et 
al, 1995]. Energetic electrons, buy virtue of their high ve- 
locities, will also be used as remote tracers of magnetic 
topology. Small X-lines should develop small-scale 
plasmoids or fluxropes. DRACO can observe whether more 
than one plasmoid/fluxrope can coexist at any one time 
and whether these plasmoids merge into more global ones 
at times. 

In this manner, magnetotail constellation measurements 
will enable us to observe the plasmoid formation process, 
determine the forces which are responsible for its ejection, 
infer plasmoid dimensions, map their magnetic topology 
and magnetic connectivity with the rest of the magnetotail 
and, finally, ascertain how plasmoids evolve with time and 
downtail distance. Specific questions that DRACO can 
answer are: How 7 much magnetospheric. mass and energy 
is carried off by plasmoids? How are multiple plasmoids 
formed and released in rapid succession and how is the 
plasma sheet replenished so quickly? How are high mag- 
netic field strength cores formed within plasmoids and 
what factors determine their magnetic topology? What are 
the causes of the rapid plasmoid expansion and downtail 
acceleration immediately following their release? What is 
the nature of the plasmoid-tail interaction and how does it 
influence neutral line motion and plasmoid ejection speed? 

1.5 Nonlinear Dynamics of the Chaotic 
Magnetotail 

We have been summarizing our current knowledge of the 
magnetotail, and the insufficiency of available observations. 
In this section, we identify the magnetotail as an example 
of a nonlinear system exhibiting chaotic behavior. In non- 
linear dynamics, chaos is defined as behavior that is ap- 
parently random and irregular, but that is actually gov- 
erned by a finite set of differential equations and there- 
fore deterministic. Simple examples include the asynchro- 
nously driven pendulum, or a ball bouncing on abase plate 
that vibrates up and down to counteract dissipation and 
maintain the bounce, or a dripping water faucet. These 
systems are accurately described by simple sets of differ- 
ential equations. The systems are precisely determinate, 


yet they appear highly irregular. The magnetotail is cer- 
tainly much more complex than these simple systems, and 
much more irregular in many ways. Nevertheless, some 
limited regularities of the magnetotail are loosely analo- 
gous to a nonlinear oscillation: for example, the recurrent 
cycle of stretching, then relaxation and plasmoid shedding 
that we identify with substorms. Magnetotail dynamics is 
of course inherently nonlinear. 

Chaotic systems have the property of limited predictabil- 
ity, as distinct from the strong predictability of the small 
amplitude, synchronously driven pendulum. The linear 
pendulum, having been disturbed, eventually settles back 
to its original routine. By contrast, when the driving am- 
plitude is sufficiently large, chaotic behavior results and 
the smallest disturbance of the pendulum switches it to a 
different but equally chaotic pattern. The pendulum's be- 
havior is completely determined by the initial condition 
and forcing function, but one cannot predict its behavior 
for any considerable period without extremely precise 
knowledge of the initial drivers. The divergence between 
two slightly differing sequences grows exponentially, the 
difference growing with each successive interval by a 
roughly constant factor, so that the smallest uncertainty 
renders predictions useless after a finite period. 

This illustrates the concept of a "horizon of predictability,” 
shared by all chaotic nonlinear systems. In practical terms, 
the horizon of predictability is an impassable barrier. At- 
mospheric movements, governed by well-defined equations, 
are thought to be in a state of chaos. Improved knowledge 
of initial and boundary conditions has significantly im- 
proved atmospheric predictability in recent decades as 
improved space-based measurements have converged with 
huge improvements in computational capabilities. Ulti- 
mately. however, incremental improvements come up 
against the exponential divergence of minutely different 
initial conditions, and the “horizon of predictability" reaches 
a practical limit. Certain statistical features of climate, such 
as annual cycles of temperature and rainfall, are exempt 
from the limitations of chaos. Other large-scale processes 
may still allow long-range prediction: the more detail one 
asks for in a forecast, the sooner it will lose its validity. 

The magnetotail has to date been studied at a primitive 
level in comparison with our lower atmosphere. Current 
observational knowledge of the magnetotail is perhaps akin 
to the meteorological knowledge of the early explorers, 
made from individual ships, widely separated in time and 
space. In contrast, theoretical work on the magnetotail is 
advanced, albeit relatively untested through confrontation 
with observations. One key motivation for the magneto- 
spheric constellation mission is to bring the tools and les- 
sons of nonlinear dynamics to the study of the magnetotail 
system. The foremost of these that is lacking at present is 
the measurement of initial and boundary conditions with 
resolution and accuracy commensurate with current and 
future generations of global simulation models. 
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1.6 Science Objectives and Questions 

The scientific foundation of the Magnetospheric Constel- 
lation DRACO mission is the quest for an understanding 
of the nonlinear dynamics of the structured magnetotail. 
This understanding will flow from answering the question: 

“How does the dynamic magnetotail control energy 
flow?” This can be broken down into derivative questions: 

• What processes control magnetotail structure and dy- 
namics ? 

• How do magnetotail dynamics respond to extrinsic con- 
ditions as well as to internal instabilities? 

• How do physical processes and regions couple over the 
hierarchy of scales? 

A more detailed expansion of these questions under mis- 
sion objective headings is given below. 

1.6.1 Determine Dynamics 

• Are there states approximating equilibria of the 
magnetotail? 

• Does the magnetotail fluctuate strongly even when the 
solar wind is steady? 

• How do these fluctuations depend on the pressure, 
speed, magnetization, and fluctuations of the solar 
wind? 

■ What is the nature of global magnetotail 
reconfigurations? 

• How do magnetotail dynamics accelerate and transport 
energetic particles? 

DRACO will make it possible to identify both equilibrium 
states and their fluctuating components as dynamic 
sequences. 

1.6.2 Understand Responses 

• Are magnetotail flow vortices generated spontaneously 
or in response to solar wind disturbances? 

• Is the occurrence of kink, ballooning, tearing or recon- 
nection instabilities determined by magnetotail bound- 
ary conditions? 

• Are fast Earthward flows and associated currents the 
same disturbances that form and release plasmoids, 
and are they reliably triggered by specific solar wind 
variations? 

• Is particle acceleration triggered or enhanced by solar 
wind discontinuity effects on magnetotail dynamics? 

DRACO will permit us to unambiguously determine and 
characterize the magnetotail's unique responses to exter- 
nal and internal drivers. 


1.6.3 Reveal Connections 

• Docs the plasma sheet collapse from the inside out 
(rarefactively), or from the outside in (compressively)? 

• How are observed ionospheric current systems and 
auroral expansions connected to plasma flows, particle 
acceleration, plasma injections? 

■ How are disturbances communicated within the 
magnetosphere? 

■ What is the horizon of magnetotail predictability at 
present and can it be extended with improved knowl- 
edge of initial and boundary conditions? 

DRACO’s extended measurement network will provide the 
first vector field movies of the magnetotail, revealing the 
connections of causes-to-effects via both convecting and 
propagating disturbances. 

1.6.4 Summary 

Magnetospheric Constellation DRACO will for the first 
time clearly identify the boundary conditions that are as- 
sociated with configuration instabilities of the magnetotail. 
In contrast with MMS, which will examine the role of local 
velocity space micro-instabilities, DRACO will provide a 
test of the macroscopic conditions for instability. 

In summary the overarching goal of the Magnetospheric 
Constellation DRACO mission is to reveal magnetotail dy- 
namics, its responses to external conditions determined by 
the solar wind, and the connections within the magnetotail 
and to its boundaries. To achieve observational closure with 
theory-based global circulation models, leading to confi- 
dence in our physical description, wc must observe the 
plasma sheet system as a coupled whole. This requires a 
network of measuring stations dense enough to resolve the 
principal vector field features of the magnetotail as they 
evolve under varying solar wind conditions. 

2.0 Science Requirements 

2.1 A New Leap in Magnetospheric 
Physics 

This section provides a derivation of Magnetospheric Con- 
stellation DRACO mission requirements from the scien- 
tific objectives, as stated in the previous sections. A con- 
stellation mission is a fundamental departure from the 
prior state of the art in space physics. In the past, space 
exploration was conducted first from lone spacecraft and 
later from small numbers of widely separated spacecraft 
(i.e., ISTP), seeking to understand the gross coupling be- 
tween principal regions. Now we have observations from 
small groups of spacecraft flying in tetrahedral formation 
(CLUSTER), examining the detailed structure and dynam- 
ics of specific regions such as the sites of reconnection 
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(MMS). Imaging missions like IMAGE, and others like 
TWINS, will give us global perspectives for emissive phe- 
nomena that can be characterized as scalar image fields. 

A constellation mission will provide a qualitative and quan- 
titative leap beyond these techniques to a distributed net- 
work of identical observing stations, providing dense vec- 
tor field sampling of a large region of space. The constella- 
tion approach flows naturally from the history of dynamic 
meteorology, where a similar approach has allowed the 
testing of comprehensive theoretical models and made re- 
liable prediction possible. It fits logically into the hierar- 
chy of physical scales in space and completes the corre- 
sponding strategic hierarchy of mission architectures. It 
implies a set of considerations for data handling, visual- 
ization and synthesis that were heretofore unique to at- 
mospheric meteorology, but now must be dealt with for 
future progress in space physics and the Sun-Earth 
Connection (SEC). 

In this section, we show how the science objectives of the 
previous section lead inevitably to the requirement for a 
constellation mission. The deployment of a constellation 
of spacecraft will be most challenging in the area of re- 
sources. In place of the conventional stress on the develop- 
ment of unique and highly capable measurement instru- 
ments, a constellation must be based on techniques more 
akin to those of the consumer electronics industry. A vari- 
ety of highly sophisticated devices must be designed into a 
durable and reliable package that can be reproduced and 
operated inexpensively The art of designing a constella- 
tion mission then clearly lies in the definition of a highly 
capable but also cost-effective instrument package and 
spacecraft, capable of achieving the science objectives of 
the mission through deployment in unprecedented num- 
bers. For any conceivable resource envelope, a trade clearly 
exists between the capability of each individual spacecraft 
and the scope and resolution of the constellation that can 
be deployed. We identify the measurement requirements 
for each spacecraft, then the required deployment density 
and domain. Finally, we identify the required instrument 
payload and spacecraft capabilities, implicitly prioritizing 
a maximal number of spacecraft. 

2.1.1 Lessons Learned from Dynamic 
Atmospheric Meteorology 

Earth s magnetosphere is a complex geophysical system 
formed by interactions between the magnetized solar wind 
and the neutral and ionized atmospheres and magnetic 
fields of the Earth. It is certainly valid and increasingly 
common to think of the magnetosphere as a part of the 
Earth’s overall weather machine \ CaJder , 1974]. "Weather” 
may be defined simply as the cosmic dynamical phenom- 
ena that are of interest to humankind by virtue of their 
occurrence within space that we inhabit. Increasingly, we 
inhabit the space around the Earth out to the limits of our 
atmosphere and beyond, into the atmosphere of our star, 
the Sun. We also invest in the placement of robotic agents 


within that space, whose survival is important to us for 
various practical reasons. Thus, magnetospheric phenom- 
ena are beginning to have relevance to our lives and are 
becoming more familiar to us in this space age. 

A single view of the global aurora suggests strongly that it 
is fundamentally a turbulent fluid flow phenomenon, as 
shown in Figure 7. The global structure of the polar au- 
rora is comparable to that of the global cloud cover of our 
planet, as exemplified in Figure 8. Clouds and aurora pro- 
vide telltale identifications of fluid parcels with distinguish- 
ing characteristics of moisture content and particle pre- 
cipitation, respectively. As parcels are transported and dis- 
torted, they reveal much about the flow patterns that ex- 
ist within their parent fluids. Large-scale eddies and vor- 
tex flows are particularly obvious in time-lapse motion pic- 
tures of these phenomena. 

The modern account of global atmospheric winds derives 
from the work of meteorologists over the past two and a 
half centuries. In the words of Nigel Calder [ 1974], "above 
all, the modern account is based on millions of measure- 
ments of conditions high over our heads, made by balloons, 
aircraft, and satellites.” In combination with the theoreti- 
cal work of Rossby, Starr, Charney, and Lorenz (some of 
which lead to a more general theory of chaos), and their 
working computer models of global circulation, a high level 
of understanding has been reached. Te rms like jet stream, 
front, depression, pressure ridge, and even baroclinic in- 
stability, have entered the vocabulary of the general pub- 
lic through daily discussion by news media weather fore- 
casters everywhere. Clearly, modern communication tech- 
nologies, most recently geosynchronous satellite imaging, 
were enabling in the synthesis of worldwide data into me- 
teorological understanding and predictive know-how. 

In the troposphere, planetary rotation sets up a balance 
between the pressure and Coriolis forces, stabilizing 
against planet-wide Hadley convection, or thermal over- 
turning of the atmosphere. However, the thermal conduc- 
tivity of the atmosphere is insufficient to support the avail- 
able energy flux. To move more heat from the equator to 
the poles, the system develops a baroclinic instability whose 
circulation permits the energy flux to be transported 
convectively. albeit in a complex 3-D pattern involving a 
mid-latitude overturning cell in addition to polar and equa- 
torial cells. The mean meridional motions (in latitude) are 
illustrated in Figure 7. Superposed on these mean mo- 
tions are the frontal and storm phenomena that make up 
mid-latitude weather. Sustained global and local monitor- 
ing of weather conditions, combined with sound theoreti- 
cal underpinnings, have permitted weather prediction to 
become usefully accurate over time scales of up to a week 
or two, though useful climate predictions continue to elude 
us at present. 

Magnetospheric meteorology is on a trajectory similar to 
that of tropospheric meteorology. Pioneering thinkers such 
as Dungey, Axford, Hines, and Nishida developed the ba- 
sic description of the overall circulation of the magneto- 
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Figure 7. This figure is an example of the global auroral structure during a significant magnetospheric disturbance. This image is 
courtesy of the Air Force Research Laboratories and was acquired in a strip wise fashion as the Defense Meteorological Spacecraft 
Program (DMSP) spacecraft moved along a relatively low Earth orbit. 


sphere. Corresponding to the Hadley cell for the magneto- 
sphere, this model of magnetospheric circulation gets the 
overall pattern correct, but overlooks much of the impor- 
tant and interesting dynamics. It was realized almost im- 
mediately [Dungey 1966; McIIwafn , 1967; see acknowledg- 
ments] that further progress on global dynamics would be 
limited with individual observing station (spacecraft) data. 
They recognized that literally “bunches” of spacecraft would 
have to be put in orbit to definitively understand the local 
and dynamic behavior of the magnetotail and its auroral 
plasmas. Since that time, considerable theoretical work has 
been accomplished and indeed, numerical models have been 
constructed of the magnetosphere and its dynamics, as ex- 
emplified in Figure 1 . Nevertheless, the small number of 
multipoint simultaneous measurements has greatly ham- 
pered closure between models and observations. The true 
behavior of the magnetosphere eludes us and we have many 
competing theoretical ideas with no way of refuting most 
of them. Consequently, we have an unrelenting debate, and 
a poor basis for understanding and predicting magneto- 
spheric dynamics. 

The magnetosphere has many storm-like departures from 
the mean expected circulation. The basic driving force, 
analogous to the thermal driving of the troposphere, is the 
solar wind flow, which seeks to impart momentum to the 
outer layers of the magnetospheric plasma. Large cycles 
of magnetotail distension and collapse suggest an insta- 
bility that might be called "barotensile" in recognition of 
the disruption of balance between plasma pressure and 
magnetic tension. Just as the true nature of the baroclinic 
instability emerged from millions of measurements of con- 
ditions throughout the troposphere, the true nature of mag- 


netospheric configuration instabilities responsible for au- 
roral substorms and storms will come from a program to 
obtain simultaneous global measurements. This will lead 
to a more complete understanding of the Sun-Earth sys- 
tem and eventually the capability to predict its behavior 
over longer time spans. We now need observations with a 
resolution in time and space that is commensurate with 
our knowledge of correlation scales in the plasma sheet, 
and of hydrodynamic wave propagation times across such 
scales: 2 R p and 10 seconds. To achieve this, we must de- 
velop affordable means by which conditions can be mea- 
sured at many points simultaneously and repetitively, and 
integrate them into space weather maps. This is the goal 
of Magnetospheric Constellation DRACO (and all other 
constellation missions). 

2.1.2 Hierarchy of Scales and Mission 
Architectures 

The progression from single satellites to constellations is 
driven by the hierarchical nature of magnetospheric phe- 
nomena. Spatial scales of magnetospheric phenomena 
range from kinetic scales associated with magnetic 
reconnection to whole-body (i.e., magnetospherewide) re- 
sponses to solar wind waves and discontinuities. Associ- 
ated with the hierarchy of spatial scales is a hierarchy of 
temporal scales, since large-scale phenomena tend to last 
longer than small-scale phenomena. Moreover, the scales 
are interrelated, one scale feeding into the next and 
vice versa. 

The hierarchical nature of magnetospheric phenomena ex- 
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Figure 8 . (Left Panel) An example of Earth’s global cloud cover imaged in visible light , illustrating the complex circulation of our 
atmosphere in response to the flow of heat from Equator to pole, the baroclinic instability. (Right Panel) Mean meridional atmo 
spheric flow, owing to the breakup of the Equator to pole cell proposed by Hadley, into three cells, leading to the formation of the Ferrel 
cell and mid latitude weather. Figures are courtesy of the University of Washington Geophysics Program. 


emplifies a general property of astrophysical plasmas. 
T. Tajima [1997] remarks that the "hierarchical nature of 
plasma manifests itself as a kind of structure made up of 
levels of physical phenomena. But these are yet ultimately 
interrelated to each other/' Tajima illustrates the point with 
the example of magnetic reconnection, a local process that 
occurs as a result of global dynamics and in turn gives rise 
to global dynamics. 

Meteorological phenomena also exhibit a hierarchy of 
scales, which feed into each other. The analogy to meteo- 
rology helps explain why magnetospheric physics must 
advance by progressing from single satellites to constella- 
tions. Meteorology has approached the problem of under- 
standing atmospheric phenomena by formulating the equa- 
tions of motion of the atmosphere and developing an ob- 
servational procedure for measuring the terms that enter 
these equations. The most important of the dynamical 
terms are wind, pressure, and temperature. For the most 
part, these quantities are invisible in the sense of being 
inaccessible through remote sensing. Meteorologists have 
found that the way to advance in the face of a hierarchy of 
dynamical scales driven by invisible forces is to establish 
an in situ observing network. The meteorological observ- 
ing network has grown to become dense enough and broad 
enough to display atmospheric phenomena on all relevant 
scales. It is also broad enough to be able to follow weather- 
producing atmospheric phenomena through their life 
cycles. To achieve an adequate network, meteorology moved 
from single isolated stations to networks in w'hich the num- 
ber of nodes grew to encompass the phenomena. 


Magnetospheric physics is in a closely analogous situation 
to meteorology with hierarchical dynamical modes and in- 
visible dynamical parameters, which in this case also in 
elude electromagnetic fields. In this type of situation, the 
problem of data-theory closure has no known solution other 
than that developed in meteorology — a dense, broad net- 
work of in situ observing stations. In the magnetospheric 
context, this means moving from single satellites to net- 
works of satellites — single satellites, then clusters, then 
constellations — until the magnetosphere itself is 
encompassed. 

Figure 9 schematically illustrates the hierarchy of mag- 
netospheric scales and how these scales relate to upcom- 
ing multispacecraft missions in the study of magneto- 
spheric phenomena (Magnetospheric Constellation 
DRACO and Magnetospheric Multiscale). Characteristic 
spatal scales on the vertical axis increase from microscale, 
mesoscale, and macroscale, to global, spanning physical 
dimensions of 10 km to 40 R f . These are deliberately vague 
to allow many phenomena to be included. "Microscale” 
roughly represents kinetic (nonmagnetohydrodynamic or 
velocity space) processes. “Mesoscale’’ represents phenom- 
ena that have scale sizes comparable to the internal struc- 
tures of the magnetosphere. "Global'’ phenomena engage 
the whole magnetosphere. Scales for duration are repre- 
sented by phenomena having those durations. Time scales 
for duration increase on the horizontal axes ranging from 
10’s of milliseconds to hours. 
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Figure 9. Hierarchy of scales in the magnetotail, illustrating the complementary 
design of various mission architectures in the study of representative phenomena. 


Figure 9 indicates key processes and phenomena that 
MMS and Magnetospheric Constellation DRACO will ad- 
dress, demonstrating the different regions of space/time 
they will resolve. All scale sizes and durations are repre- 
sented and phenomena at different scale sizes are related 
by mode of origin. This illustrates for the magnetospheric 
case the point that Tajima made about astrophysical 
plasma in general: They are hierarchical and their scales 
relate to each other. To understand the system fully, we 
must resolve all size/time scales systematically and, ide- 
ally simultaneously. 

At present, single satellites have only glimpsed the phe- 
nomena shown in Figure 9, and have not been able to 
resolve inherent space-time ambiguities. As such, single 
satellites only describe the state of phenomena that change 
slowly, such as the relatively slow-moving magnetospheric 
boundary layers, inverted V’s, and diffusive changes in the 
ring current and radiation belts. For other phenomena, 
single satellites only give statistical information accumu- 
lated by multiple encounters with the phenomena on sepa- 
rate occasions. In meteorological terms, for other phenom- 
ena single satellites give climatologies rather than descrip- 
tions of actual weather systems. 

CLUSTER-class missions such as MMS will be able for 
the first time to resolve the state and evolution of the 


plasma physics of reconnection and turbulence down to 
microscopic scales (10 km). At the other end of the MMS 
range, it will also capture appreciable elements of such 
phenomena as substorm triggers and bursty flows. But 
there are important science questions that lie outside the 
portion of parameter space that presently planned clus- 
ter-class missions can explorer. 

As outlined in the previous sections, many important pro- 
cesses are operating in the constellation-class domain cov- 
ered by Magnetospheric Constellation DRACO. These in- 
clude the magnetosphere s best-known and most energetic 
dynamical modes: the substorm. In addition, for instance, 
it is within this constellation-class domain that the global 
consequences of reconnection can be understood. Other con- 
stellation-class processes include manifestations of flux 
ropes. 

In summary, the hierarchical, interdependent phenomena 
of astrophysical plasmas in general and of magnetospheric 
plasmas in particular drive a progression from single sat- 
ellites to constellations. Meteorology has similar dynam- 
ics and has implemented a corresponding progression of 
atmospheric observing stations. Magnetospheric physics 
has passed through a single-satellite state and is on the 
brink of a cluster state. The time is right to plan for the 
constellation stage. 
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2.1.3 Complementary with Other Missions 

NASA SEC missions are coordinated to provide knowledge 
about the magnetosphere at all relevant scales. The ISTP 
missions provided the most global coverage to date, estab- 
lishing the importance of continuous solar and solar wind 
monitoring in conjunction with single -spacecraft observa- 
tions at other strategically chosen locations. In the pro- 
cess, ISTP has also demonstrated the capability to remotely 
sense and image the energetic ions of the inner magneto- 
sphere [Henderson etal., 1997, Jorgensen etal, 1997]. Now, 
the IMAGE and TWINS missions are putting in place the 
tools for routine dedicated imaging of the inner magneto- 
spheric energetic ions as well as the cold plasmaspheric 
helium, low-energy ionospheric outflows, and the motions 
of the dayside magnetopause. However, imaging missions 
are unable to probe the magnetotail region, owing to the 
low plasma and neutral gas densities in that region. Even 
with increased sensitivities, magnetospheric imaging mis- 
sions cannot make vector measurements of the magnetic 
field and plasma flow field, being limited to scalar quanti- 
ties such as line-of-sight flux. These are extremely impor- 
tant in the magnetotail region. 

The MMS mission will focus upon the microscale and ve- 
locity space processes that occur in magnetospheric bound- 
ary layers, including the thin magnetotail current sheet: 
reconnection, particle acceleration, and turbulence. These 
processes regulate the flow of mass, energy, and momen- 
tum from magnetotail lobe storage regions into the inner 
magnetosphere and ionosphere via high-speed plasma 
flows, energized particles, and field-aligned currents* 
thereby determining when and where geomagnetic distur- 
bances will occur. MMS will discriminate between spatial 
and temporal phenomena and provide the local gradient 
observations needed to characterize and then identify these 
poorly understood microscale processes. Fundamentally, it 
will address the influence of microscale velocity space struc- 
ture on the dynamics of the magnetosphere. 

The Geospace Electrodynamics Connection (GEC) mission 
will define the nature of the connections linking the mag- 
netosphere and the ionosphere. The key objectives of this 
mission include: 

• Identifying the spatial and temporal scales over which the 
ionosphere and thermosphere respond dynamically to mag- 
netospheric input, 

• Establishing the route linking magnetospheric and iono- 
spheric currents as a function of conductivity, and 

• Defining how plasma flows between ionosphere and mag- 
netosphere affect the dynamics of both regions. 

Whereas the MMS and GEC missions will provide a wealth 
of information concerning the small-scale structure of pro- 
cesses in magnetospheric boundary layers and the iono- 
sphere, neither mission will explore the global equilibrium, 
response to disturbances, and macroscopic configuration 
space instability of the magnetotail, nor its global link- 


ages. Specifically, DRACO will focus on the global conse- 
quences of the reconnection process, the global control of 
its onset locations and trigger mechanism and evolution, 
and the causal links of (possibly) multiple reconnection 
sites. Only the dense network of vector measurement sta- 
tions provided by DRACO can address these fundamental 
global issues, complementing the other STP missions. Fun- 
damentally, DRACO will address the influence of macroscale 
configuration space conditions on the dynamics of the mag- 
netosphere, complementing the other missions. All of the 
missions are described in web pages at http:// 
s t p . gs fc. n asa . gov 

2.1.4 Data Synthesis Considerations for 
Constellation-Class Missions 

2. 1.4.1 Novel Modeling, Data Handling and Assimi- 
lation Techniques. The Magnetospheric Constellation 
DRACO mission will provide multipoint in situ magneto- 
spheric measurements of substantially different quality 
and quantity than any previous mission. The quality of 
the measurements lies not so much in the quality of the 
individual measurements — which are indeed comparable 
to previous missions — but in the fact that these observa- 
tions are made in a coherent fashion on a large number of 
spacecraft simultaneously. Making optimal use of the mea- 
surements will require the development and application 
of new techniques. The implementation of the DRACO mis- 
sion will result in a leap, similar in its importance to that 
made by the atmospheric physicists many decades ago. In- 
dividual weather observations were gathered into synop- 
tic maps (summarizing a large number of individual ob- 
servations). Later, the maps came to be further processed 
using spatial objective analysis and data assimilation. The 
analysis of individual spacecraft observations will become 
mostly obsolete and will be replaced by synoptic analyses 
that use observations from all satellites in the constella- 
tion. An example of the type of information that will be 
obtained from DRACO is shown in Figure 10, generated 
from the global simulation and visualization work of 
Goodrich et al., [1999]. 

Figure 10 shows a snapshot of the magnetotail neutral 
sheet region from a global simulation, showing the forma- 
tion of a strong Earthward flow in the evening sector. Ar- 
rows indicate the local flow vector, while color contouring 
indicates the flow of electric current into or out of the neu- 
tral sheet region, along magnetic field lines. 

Techniques developed for the atmospheric sciences over 
the past decades can be borrowed to solve the new analy- 
sis challenges. This is critical for the success of DRACO 
because — as the atmospheric sciences show us — the de- 
velopment of these techniques can take decades and some- 
times lead into dead ends. By exploiting the existing knowl- 
edge, DRACO will make these tools ready for the first 
magnetospheric multipoint mission in less than a decade. 
It is crucial for DRACO s success to include a strong mod- 
eling element right from the beginning. These analysis and 
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Figure 10. A single spacecraft snapshot gives good information at a point in space and time. A single spacecraft orbit gives good 
information but spread over a 3-day orbit. A constellation gives a dense set of data with good time resolution and spatial resolution 
adequate to resolve the main global features of the simulation fright panel). Figure is courtesy ofC. C. Goodrich , from Goodrich et al. 
[ 1999j . 


assimilation tools need to be ready when the first DRACO 
measurements become available. The development of these 
tools will require substantial effort. 

The first factor is the simultaneous data coverage of the 
magnetosphere. A significantly nonuniform distribution of 
the satellites will result in a loss of information in regions 
with sparse coverage, while an excessively dense concen- 
tration of the space probes in relatively small regions will 
not enhance the fidelity because of limited spatial resolu- 
tion of the mathematical models. Therefore, the inversion 
methods must include a means of "filling in the gaps”; for 
example, using the predictions of the global statistical 
models in regions with poor spatial coverage. Similar prin- 
ciples have been successfully applied in consolidating 
multipoint ionospheric measurements in the Assimilative 
Mapping of Ionospheric Electrodynamics (AMIE) models. 
In this regard, simultaneous monitoring of the conditions 
of the incoming solar wind could greatly enhance the ac- 
curacy of the DRACO-based modeling. In spite of the fact 
that the DRACO measurements will be made mainly at 
low latitudes, the dynamic evolution of the orbits and dis- 
tortions of the tail will provide sufficient information about 
structure normal to the orbital plane [Tsyganenko, 1998]. 

The second factor is the flexibility of the model represen- 
tation; that is, its ability to reproduce a large variety of 
possible magnetic configurations. Substantial progress in 
this direction has taken place recently. A mathematically 
simple and general approach was developed to model the 
largest-scale geomagnetic field, produced by currents on 
the magnetospheric boundary. This magnetopause contri- 
bution is important throughout the magnetosphere and 
depends mainly on the dynamical pressure of the solar 
wind. The cross-tail electric current, flowing from dawn to 
dusk in the equatorial plasma sheet, is also critically im- 
portant. The most recent mathematical models [e.g., 
Tsyganenko and Peredo, 1994], superpose the fields from 


several equatorial current sheets with variable thickness 
and different radial profiles, all confined within a magne- 
topause with realistic shape and size. In addition, power- 
ful new techniques became available recently making it 
possible to take into account the seasonal/diurnal warping 
of the tail current and its twisting around the Sun-Earth 
line, induced by the reconnection of the geomagnetic and 
interplanetary magnetic fields [ Tsyganenko, 1 998]. Now 
many more degrees of freedom can easily be introduced 
into the model, making it possible to magnetically map 
the effects of catastrophic reconfigurations of the 
magnetotail during disturbed periods. Further progress in 
this area is expected on the basis of the new deformation 
method \ Tsyganenko, 1998], allowing one to greatly expand 
the variety of modeled magnetic configurations, including 
plasmoids and other transient features. 

2. 1.4.2 Data Assimilation in Dynamic MHD 
Modeling. Data assimilation takes the synoptic mapping 
a significant step further by introducing the "first prin- 
ciple” dynamical MHD models into the analysis process. 
The underlying idea is that assimilation of data can mini- 
mize errors in both the data and the model, even when the 
model suffers from simplifying assumptions, and thus pro- 
vide an optimal "map” of the physical state of the mag- 
netosphere at any given time. 

Data assimilation techniques have been successfully used 
in atmospheric and oceanographic models for some time 
\Ghil et ah, 1997], but never before in magnetospheric 
modeling. Data assimilation techniques take the observa- 
tions distributed in time and space and merge them dy- 
namically into the numerical model in order to determine 
as accurately as possible the instantaneous state of the 
magnetosphere. The primary objectives and needs for the 
data assimilation are our poor knowledge of the initial state 
that is used to start a model simulation and the fact that 
the simulated magnetosphere may gradually diverge in 


18 


Magnetospheric Constellation Mission... 


time from its actual state. Data assimilation ensures that, 
in the course of time, the state of the simulation remains 
optimally close (in terms of a merit function) to the con- 
current observations at the measurement points. The pri- 
mary reason that data assimilation methods have not yet 
found their way into magnetospheric modeling is the scar- 
city of magnetospheric in situ measurements. The Mag- 
netospheric Constellation DRACO mission will change that 
dramatically by providing for the first time a fairly dense 
grid of in situ data describing the MHD state of the entire 
modeled system. The fact that these observations are re- 
stricted to the vicinity of the tail plasma sheet is no major 
obstacle. Atmospheric models have to deal with similar 
unevenly distributed data sets and powerful methods ex- 
ist to make optimal use of them [Miller et al, 1994; 
Talagrand, 1997]. 

The power of data assimilation has been demonstrated con- 
vincingly by the numerical weather prediction models which 
all use assimilation techniques. Employing these techniques 
in magnetospheric modeling will result in a dramatic break- 
through, since the data assimilation will serve as an ulti- 
mate test of the underlying model’s quality. A systematic 
divergence of a model from the data will provide a clue to 
some intrinsic deficiency to be corrected. In addition, not only 
will the model improve overall, but also new knowledge will 
be gained about the physical processes, whose incomplete or 
wrong description had caused the model to fail. In short, the 
data assimilation forces a continuing reevaluation of our 
understanding of how the magnetosphere w r orks. This is one 
of the primary reasons why the DRACO Mission should be 
implemented. 

2.2 Mission Requirements 

2.2.1 Measurement Requirements 

The science objectives for DRACO are expressed in terms 
that identify the most fundamental measurements: mag- 
netic field, plasma flow field, and energetic particle accel- 
eration. These are the three basic measurements that 
DRACO will require. It may be objected that charged par- 
ticles are accelerated by electric fields, so that electric field 
measurements are required as well. However, at the time 
scales of interest to DRACO (~ 10 sec), and in the region of 
interest to DRACO, the electric field may adequately be 
inferred from the plasma flow (transverse component) and 
from measurements of magnetic field variations (induced 
parallel component). Short time-scale phenomena can usu- 
ally be best observed using electric field measurements to 
integrate instantaneously over the plasma velocity distri- 
bution. However, such phenomena are generally beyond 
the scope of DRACO and can be best diagnosed by mis- 
sions devoting their resources to detailed microphysical 
observations. Table 1 provides a compact summary of the 
measurement requirements for DRACO. A translation of 
these requirements into a straw man payload is provided 
in this section. It is the nature of DRACO science objec- 
tives to require coordinated measurements of all MHD 


plasma state variables. These individual measurements 
cannot be exclusively associated with specific DRACO ob- 
jectives; all contribute equally. 

2.2.2 Orbit-Attitude Requirements 

2.2.2. 1 Orbit. DRACO must sample the magnetotail re- 
gion from about 7 R E to 40 R £ in the night side along the 
Sun-Earth line, and across the magnetotail east and west 
of the Sun-Earth line by - 10 R £ . This is based on the im- 
portance of phenomena within the region from geosynchro- 
nous orbit to beyond the convection reversal associated with 
substorm activity, near 25-30 R £ . Eccentric orbits with peri- 
gees near 2-3 R £ provide an optimal trade between orbit 
stability and economical transmission of commands and 
data. A set of orbits with apogees ranging from 7 to 40 R £ 
would provide all the coverage that is needed. The number 
of spacecraft required must be determined on the basis of 
detailed orbit modeling, which is reported below in the 
mission design section. 

2. 2. 2. 2 Attitude. Spacecraft attitude stability and con- 
trol are required to within 5°. The attitude stability re- 
quirement is derived from the mission requirement for 10% 
absolute knowledge of the magnetic field. Data collection 
in shadows and shortly thereafter is not required, but it is 
highly desirable. Since thermal changes during shadows will 
affect the magnetometer and particle measurements, as spin 
rate changes, it is recommended (but not required) that spin 
rate data be obtained interchangeably from the Sun sensors 
and the Earth/Moon sensors (infrared detector). 

2.2.3 Spatial Resolution and Number of 
Spacecraft 

DRACO must provide a typical interspacecraft spacing 
(resolution) of ~ 2 R £ . Moreover, the distribution of inter- 
spacecraft spacing should be as narrow as practical, so that 
random "holes" and "bunches” in the spacecraft distribu- 
tion are minimized. This requirement is translated into 
orbital strategy in Section 3. 

2.2.4 Mission Duration 

To achieve its science objectives, DRACO must have a mis- 
sion lifetime of sufficient duration to observe magnetotail 
behavior during a variety of solar wind conditions and at 
least a sampling of typical solar wind disturbances, includ- 
ing co-rotating interaction regions and Coronal Mass Ejec- 
tions (CME's). The entire constellation should be put into 
orbit with apogees in the local time range around local 
dawn. Initial activation would then be complete well in 
advance of the rotation of the tail across the constellation 
orbit pattern as the Earth orbits the Sun. This would pro- 
duce a pioneering but ancillary study of the dawn flank of 
the magnetosphere and its interaction with the 
magnetosheath, early in the mission. This would be fol- 
lowed by a main mission phase as the constellation passes 
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Measuremen t 

Range 

Resolution 

Time Resolution 

Comments 

3 Axis 

Magnetic Field 

+/-300nT 

0.1 nT 

1 sec 

I min @ 100 Hz 

Fluxgale technology exists 

Piasma 2D Tern 
perature 

10- 2000 eV 

20% 

10 sec 

Requires 18ff field of view 

Plasma Flux 

Iff - iff 
cm 2 s l sr 1 
eV/eV 

20% 

10 sec 

Electrostatic analyzer technology' exists 

Plasma 3-D Velocity 
Electron PAD 

i - WOO km/ s 

20% 

2ff 

10 sec 

Mass analysis significant but not absolutely 
required 

Particle Energy 

20 - 500 keV 

20% 

1 0 sec 

Mass analysis significant but not absolutely 
required 

Particle Flux 

IxlCP 1x1 ff 
cm 2 s' 1 sr* 

20% 

10 sec 

Solid-state telescope technology r exists 

Particle Pitch Angle 

ISff 

2ff 

10 sec 



Table 1 . Measurement requirements : 


through the center of the magnetotail over a period of 
roughly 4 months. This, in turn, would be followed by a 
study of the dusk flank and boundary layer interactions. 
All this would occur during the first 6 months of opera- 
tion, and would achieve the minimum success criteria for 
the mission. Continued operations over a subsequent year 
would produce unprecedented studies of the day side mag- 
netopause, magnetosheath, bow shock, and foreshock re- 
gions in the solar wind. Studies of these regions are of 
course considered to be ancillary science for DRACO, but 
would have high value, and would position DRACO for a 
second pass through the tail 1-2 yrs after launch. 

2.2.5 Data Volume and Flow 

The general data flow scenario calls for data from an en- 
tire orbit to be stored on board, and then downloaded 
quickly during perigee passes near the Earth. This will 
minimize the requirements for transmitter power and 
ground operations, but creates the need for onboard stor- 
age. Table 2 presents the Magnetospheric Constellation 
DRACO science data accumulation rates and total data 
storage requirements per orbit, absent any data compres- 
sion. Loss-less compression will be used routinely in prac- 
tice, creating a factor of two contingency in the overall data 
volumes shown here, reducing the requirement for bulk 
memory on the spacecraft, and possibly reducing the down- 
link requirement from that reported here. The typical in- 
trinsic instrument sampling rates will generate the acqui- 
sition of a quantity of data from a single spacecraft that is 
quite manageable by contemporary standards. 

Data assimilation is greatly facilitated by homogeneity of 


the individual data sets in terms of time, energy, and an- 
gular resolution. Thus, the economic advantage of making 
all the spacecraft and instruments identical is synergistic 
with the overarching science goals for the mission. In place 
of time-variable burst mode data-taking, routine data prod- 
ucts will include both high time resolution integrated data 
(moments), and lower time resolution detailed data. 

DRACO's lower apogee orbits have reduced storage require- 
ments in accord with the shorter orbital periods, the low- 
est apogee orbit having a period of ~ 1 2 hrs. Assuming iden- 
tical spacecraft, the excess data storage capability in the 
lower orbits translates into considerable flexibility in the 
planning of downlinks. All of the lower apogee data can be 
transmitted during considerably shorter contact periods, 
or it can be accumulated for multiple orbits and transmit- 
ted using the same contact period as the higher apogee 
orbits. The opportunity also exists for triggered burst mode 
data to be collected at higher than normal time resolution, 
with the burst mode data stored in the contingency memory 
for the lower apogee orbits. Provision for burst mode col- 
lection can also be made for the higher apogee orbits when 
data compression is implemented, particularly if Tossy” 
compression of detailed data is used, as is becoming in- 
creasingly accepted for image data. 

A summary of a possible scenario including burst data col- 
lection is shown in Table 3. The lowest apogee orbits can 
easily collect an amount of burst data per orbit equal to 
the routine data collection, assuming a data dump can take 
place once per perigee pass. Lower tracking availability 
will reduce the burst collection from those orbits. The high- 
est apogee orbits (periods greater than 50 hours, i.e., apo- 
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Instrument 


Basis (3 sec spin period ) 


Data Rate 


Magnetometer 

(MAG) 

16 vectors/spin @ 16bits/sample 

768 bits/spin 

Plasma velocity ; 
ions 

8 elevation x 11 azimuth x 16 energy- x 8 bits/sample 

11264 bits/spin 

Plasma velocity ; 
electrons 

8 elevation x 11 azimuth x 16 energy x 8 bits/sample 

1 1264 bits/spin 

Energetic ions 

4 elevation x 8 azimuth x 16 energy x 8 bits/sample 

4096 bits/spin 

Energetic electrons 

4 elevation x 8 azimuth x 16 energy x 8 bits/sample 

4096 hits/spin 

Total/Spin 


31488 bits/spin 

Total / 1 00 hrs. 

Maximum apogee orbit (40 RJ data accumulation 

3.8 Gbit/orbit maximum 


Table 2. Instrument data sampling rates. 


gees greater than 20 R E ) must retain most of the memory 
available for routine data collection. However, they can eas- 
ily accommodate an amount of burst data similar to that for 
the lower apogee orbits, simply by foregoing routine collec- 
tion in the lower altitude parts of the orbit (r < 7 R E ). 

2.2.6 Command Uplink 

Instrument commanding will be minimal. All commands 
are time-tagged and uploaded to an onboard buffer in the 
Command and Data Process Unit (CDPU) for delayed 
implementation. Commands include clock updates, provid- 
ing information on when to start and stop routine data 
collection, when to anticipate ground contact next (based 


on updated contact times), what criteria are to be used to 
trigger burst mode, and which quantity to burst on. Other 
options may include changes in the number of burst inter- 
vals, software uploads, and onboard processing options, but 
do not constitute standard commands during routine track- 
ing sessions. In general, instrument mode changes will not 
be part of the routine commanding unless required for 
health and safety. 

During initial activation, high-voltage instruments will be 
ramped and checked out during perigee pass contact peri- 
ods, then turned over to automated operation for the rest 
of their lifetimes. While the first few spacecraft activations 
may be watched closely this operation will become quite 
routine thereafter. 


Orbit 

Apogee = 7 R e 

Apogee - 40 R £ 

Orbital Period (OP) - 

15 hrs. min. 

120 hrs. max. 

Routine Accumulation 
Interval (RAT) 

12 hrs. min. 

1 00 hrs. max. 

Burst Accumulation Interval 
(BAI) 

12 hrs. (equiv. routine data) 

20 hrs. (equiv. routine data) 

Total collection per orbit 

128 MByte 

640 MByte 


Table 3. DRACO Downlink Summary 
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2.3 Instrumentation Requirements and 
Strawman Payload Resources 

The science requirements can be fulfilled with three prin- 
cipal measurement instruments: 1) Triaxial MAG 2) Elec- 
trostatic Analyzer (ESA) and 3) Solid State Telescope (SST). 
The instruments are based on current technology. Antici- 
pated advances in technology can reduce the instruments’ 
size and weight; such advances will increase the number 
of probes that can be released on a single launch, increase 
the capability for the same weight and power, and enhance 
the science return. 

The unique aspects of the science payload for DRACO are: 

• The high level of integration of instruments; 

• The high level of integration of the spacecraft around 
the science instruments; and 

• Special design considerations that will render the in- 
strument suite manufacturable in large numbers and 
capable of rapid testing and calibration as a single unit. 

The following discussion is a summary of instrument re- 
source requirements. These are derived from current tech- 
nology, with the addition of a strategy for integration to 
reduce weight and power down to a minimal amount. The 
strategy is to implement most of the electronics design on 
a single board of the spacecraft Data Processing Unit 
(DPU). Table 4 summarizes those resources for the three 
instrument types. 

2.3.1 Magnetometer 

Triaxial Huxgate magnetometer capabilities that satisfy 
science requirements have been demonstrated in previous 
missions (Lunar prospector, FAST, WIND). Careful design 
for the DRACO mission can also reduce weight and power 
to make the magnetometer commensurate with 
nanosatellites. Ring cores of a few centimeters have been 
tested in the laboratory and low-weight copper windings 
have been flown previously (POL AR/Magne tic Field Instru- 


ment (MFE)). A 16-bit digitization is sufficient. Noise lev- 
els based on previous flight experience and lab tests are 
+0.05 nT. Two ranges are necessary: 1) a low field range 
(<1000 nT) for science and is noise-limited, and 2) a high 
field range (< 1 6000nT), for backup attitude knowledge and 
ground testing of the unit but where digitization is limited 
(0.05 nT). 

The mission will benefit from a Quality Assurance (QA) 
program that emphasizes the redundancy built into the 
large number of nanosatellites rather than in a single in- 
strument realization. This makes instrument-command- 
ing and mode-switching simpler and the instrument lighter. 
In accordance with this philosophy, a single drive (rather 
than the traditional triple-redundant) electronics circuitry 
can be used, which further reduces the power and weight 
of the instrument. 

Early resolution of magnetic cleanliness issues has proven 
to be commensurate with a low-cost, high-performance 
flight unit. All known sources of magnetic material (in- 
cluding the possibility of broom magnets on the SST) must 
be compensated for early in the program. A clear under- 
standing that there is no expectation for useful data col- 
lection during Radio Frequency (RF) operation, must be 
addressed. A primary consideration for magnetic cleanli- 
ness is the current from solar panels which demands care- 
ful electrical design and wiring. 

Spacecraft shadow operations impart a change in the spin 
rate due to thermal contraction of the spacecraft body. 
Shadow data are not mission-critical, but occur in regions 
where the reconnection or current disruption processes 
erupt. Salvaging such data should be given careful consid- 
eration. One possibility is to use an Earth Infrared (IR) 
sensor to accurately track the spin rate, in addition to the 
Sun sensor. Automatic switching to the IR sensor once the 
Sun sensor signal degrades (shadow) will accomplish this 
goal at little added expense (or at no added expense if the 
IR sensor is already part of the Attitude Control Subsystem 
(ACS) system). 


Instrument 

Type 

Dimensions 

(mm) 


Mass 

(g) 

Power 

(mW) 

Comments 

MAC 

Sensor 10x10x10 
Boom: 1000x20 
Electronics : 76x96x40 


200 

300 

Sensor, High Voltage Power 
Supply (HVPS) and analog 
electronics 

Plasma Velocity 
Analyzer (PVA) 

Sensor: 100x100x100 
Electronics: 50x 1 00x50 

+90x + 10 
+90x + 10 
from radial 

1200 


Sensor, boom and analog elec- 
tronics 

Energetic Particles 
Analyzer (EPA) 

Sensor: 100x100x100 
Electronics: 50x100x50 

+90x+W 
from radial 

800 


Sensor, HVPS and analog elec 
ironies 

CDPU 

Logic: 90x90x50 
Power: 90x90x50 

N/A 

500 

600 

Digital logic and power con- 
trol 


Table 4. DRACO instrument resource allocations. 
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Figure 11. Lunar Prospector magnetometer boom latch. 


A small (30 cm) boom deployable via centrifugal force 
shortly after release is expected to be necessary for clean- 
liness purposes. The weight of the boom (from carbon ep- 
oxy based on Lunar Prospector, FAST and POLAR heri- 
tage) is 100 grams and has been included in the magne- 
tometer weight given above. An arm boom with a latch 
would be preferable to a telescopic boom, since it better 
maintains the MAG mounting orientation after release. 

Figures 11 and 12 show the Lunar Prospector boom latch, 
and the ROSETTA lander magnetometer respectively. The 
DRACO sensor can be - 1 .5 cm or smaller. Developing tech- 
nologies will make possible even smaller and lighter mag- 
netometer systems, but this is not essential to DRACO. 

Onboard calibration is a significant component of the mis- 
sion success. It is important to ascertain that two space- 
craft provide identical values when no currents are present. 
Existing methods from cross-calibration of CLUSTER mag- 
netometers can be used. These methods call for initial in- 
dividual calibration of zero levels, gains and sensor orien- 
tations, followed by inter-calibration in low-current regions 
of the magnetosphere. The benefits from such a method 
can surpass any preflight ground testing with significant 
value added to the data. Automating this procedure is es- 
sential to limit costs. 



Figure 12. ROSETTA lander magnetometer at 50g. 


2.3.2 Electro-Static Analyzer Plasma 
Instrument 

The ESA instrument measures total ions and electrons at 
energies between 10 eV and 30 keV. The science require- 
ments of Magnetospheric Constellation DRACO can be met 
with a single analyzer pair (one for ions and one for elec- 
trons) on each nanosatellite. These will collect measure- 
ments of the complete 3-D ion and electron velocity distri- 
butions during each spacecraft spin. Particles enter the 
analyzer over a - 180 degree Field Of View (FOV) and are 
selected in energy by the potential applied between the 
outer (0 Volts) and inner (0-3 kV, sweeping) concentric 
spheres. The particles are focused onto the Micro Channel 
Plate electron rmdtiplier (MCP). 

The - 1 80 degree field of view is aligned with the spin axis 
so that a full 4n steradian solid angle can be covered each 
spin. The FAST plasma instrument illustrates the flight 
heritage of hardware that would be practical to fly on 
DRACO, while the Medusa plasma instrument for the 
Swedish Munin spacecraft show r s that even smaller instru- 
mentation is being developed currently (Munin total space- 
craft mass = 6 kg). Both are based on the dual opposing 
top-hat concept. Figures 13 and 14 illustrate the assem- 
bly of 4 such ESA's for the FAST spacecraft, and the com- 
plete Munin spacecraft, weighing only 6 kg, and contain- 
ing one Medusa ESA. 

Incremental gains in weight and pow r er reduction are ex- 
pected from development funding using the top-hat instru- 
ment principle, or other related optics principles (e.g.. the 
Munin Medusa instrument). However, breakthrough in- 
strument advances may result from the development of 



Figure 13. The FAST ESA mechanical structure, 
including four independent fast plasma analyzers. 
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Figure 14. The Munin structure with solar arrays and 
a single Medusa analyzer, with a total spacecraft mass 
of 6 kg. 


delta-doped detectors and other similarly motivated solid- 
state devices. That principle is expected to reduce the 
threshold of particle detection in solid-state detectors, 
which will allow solid-state telescopes to operate in the 
plasma energy range. This is a technology that should be 
pursued further for DRACO as instrument development 
funding becomes available. 

2.3.3 Solid State Telescope for Energetic 
Particles 

The principle of solid-state telescope detection of total ion 
measurements has a long heritage in previous missions 
and can be readily implemented in a way that meets the 
science requirements of DRACO at a weight and power 
commensurate with a 10-kg nanosatellite. Using WIND 
and Equator-S heritage detectors, each double-ended tele- 
scope would consist of a set of 3 semiconductor detectors. 
The detectors are fully depleted, ion-implanted silicon with 
low leakage current. One side detector measures ions, and 
the other side detector measures electrons. High-energy 
(400 keV- IMeV) electrons are measured by the center de- 
tector. The operational principle is illustrated in Figure 
15. Ions are stopped on one side using a Lexan foil and 
electrons are stopped on the other side using a small 
"broom” magnet. A pair of such detectors stacked together 
and using a common magnet yoke is also shown in Figure 
16 (one pair is shown as deployed on WIND). Analog elec- 
tronics for pulse shaping amplifier A/D converter and co- 
incidence logic are housed close to the detector, while digi- 
tal electronics are housed at the spacecraft DPU. 



Figure 15. Solid State Telescope (SST) operational 
concept. 


Increased capability can be achieved from new technolo- 
gies utilizing miniature detector design and Application 
Specific Integrated Circuit (ASIC) technologies for analog 
and digital processing. Such capability entails full angu- 
lar resolution and mass discrimination at no additional 
weight and power expense. These ASIC units have already 
been developed and tested and are slated for flight on up- 
coming NASA missions. 


2.3.4 Signal/Noise Requirements 


Figure 17 indicates the basis for sensitivity requirements, 
using data from recent missions (Geotail and WIND) . These 
observations were obtained in the magnetotail from 8 to 
30 R e , spanning variations in flux by an order of magni- 
tude for these representative examples. The pixel geomet- 
ric factors for successful observations in this region range 
from 0.0001 cm2-sr-eV/eV to 0.01 cm2-sr-eV/eV per pixel 



Figure 16. WIND spacecraft realization of SST. 
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WIND@7R E (96/10/04) 

Geotail@9R £ (98/07/22)-— : : ; 
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WIND(g)7R E (96/10/04) — 
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Figure 17 Peak fluxes (left side scales) and count rates (inner left and right side scales, also logarithmic) are shown both from the 
LEP, or low energy' plasma instrument, and EPIC/ICS, the energetic particle instrument, on Geotail. Similar results from WIND, 
closer to the Earth, are derived from the PESA/EESA (plasma instrument), and the SST (energetic particles instrument). Both data 
points and representative fits are shown. Data is courtesy of T. Mukai (LEP), D. Williams (EPIC), and R. P. Lin (PES A/EES A/SST) 


for the plasma velocity analyzer, and from 0.0 1 to 0. 1 cm2- 
sr-eV/eV for the energetic particles instruments. 

To assure the required time resolution of 10 sec (~ 3 spins 
as currently planned) for plasma measurements having 
an accumulation time at full angular and energy resolu- 
tion of 0.06 sec (16 energies, 11 azimuths) (elevation im- 
aged on 8 anodes), we need = 18 s° count rate to have mean- 
ingful statistics. This is readily attained for the range of 
GF’s used on Geotail and WIND, as seen in Figure 17. 

To assure the required time resolution for energetic par- 
ticle measurements having an accumulation time at Full 
angular and energy resolution of 0.08 sec (16 energies, 8 
azimuths) (elevation imaged on 4 detectors), we need = 13 
s 1 count rate to have meaningful statistics. This is easily 
attained for electrons with the range of Geometric Factors 
(GFs) used on Geotail and WIND, as seen in Figure 17. 
For energetic ions, this is satisfied in the range up to a few 


hundred keV, but some integration will be required For 
the high-energy ions, with GF's similar to those on Geotail 
and WIND. 

Problematic noise sources for these measurements include 
ultraviolet (UV) leaks through the optics, and penetrat- 
ing MeV electrons. The former can be reduced adequately 
with careful blackening and baffling. The latter are only a 
serious problem inside of 10 R E , but adequate shielding 
close to the detectors must be included in the instrument 
designs, since the nanosatellite structure offers limited 
shielding external to the instrument. 

3.0 Mission Design 

The Magnetospheric Constellation DRACO mission can 
currently be implemented with 50 state-of-the-art satel- 
lites similar to NASA’s Space Tec.hnology-5 (ST- 5) space- 


Magnetospheric Constellation Mission... 


25 


craft to provide the science floor mission. NASA invested 
considerable resources in nanosatellite technology devel- 
opment, as has private industry. The first nanosatellite 
launches occur in 2000. It is expected that judicious in- 
vestments to close the gaps between current nanosatellite 
technical capabilities and the DRACO requirements will 
enable a maximum science mission with 100 satellites. 

The Magnetospheric Constellation DRACO launches from 
Kennedy Space Center on a Delta 7925 H into a 1 85 km by 
20 R e transfer orbit. A specially designed dispenser ship, 
which carries the 50 to 100 spacecraft, provides the pro- 
pulsion to raise the orbit perigee from 185 km to 3 R e and 
delivers the spacecraft to orbits with a common perigee of 
3 R e and apogees ranging From 7 to 40 R £ . 

3.1 Constellation Orbit Strategy 

The Magnetospheric Constellation DRACO must sample 
the magnetotail region from about 7 R e to 40 R E , centered 
at midnight and extending east and west of the Earth-Sun 
line by about 10 R £ . The spacing of the spacecraft should 
be as uniform as possible, given the eventual random phas- 
ing of the spacecraft around their orbits. A nested group of 
orbits with a common perigee and apogees from 7 to 40 R £ 
provides this spacing. Thus, all of the spacecraft pass 
through a common region for economical telemetry down- 
links. The 3 R £ perigee distance was chosen mainly for or- 
bital stability and communication. Careful study of the con- 
stellation orbit options has been performed to determine 
the number of nanosatellites that will give adequate spa- 
tial coverage, as outlined in this section. 

3.1.1 Orbit Studies 

As a basis for quantitatively assessing this strategy, orbit 
modeling has been performed on a suite of up to 100 vir- 
tual spacecraft. In particular, this simulation addressed 
the spatial separations obtained using various orbital de- 
ployment strategies, with an eye toward meeting the mis- 
sion measurement requirements with the minimum pos- 
sible number of spacecraft. It was determined that the 
optimal deployment strategy is a distribution of single 
spacecraft in closely spaced orbits, with random placement 
of spacecraft in orbital phase. 

Figure 18 illustrates the optimal strategy using 2-D plots 
of the probability density of the spacecraft in lxl R E cells. 
The optimal approach is compared with another strategy 
that was considered, in which families of spacecraft are 
placed in common orbits. A much more uniform distribu- 
tion results from the uniformly distributed approach. Of 
course, any orbit chosen will be susceptible to perturba- 
tions and evolve over time. 

When the constellation is examined statistically in terms 
of the distance from each spacecraft to its nearest neigh- 
bor, the result is a histogram of nearest neighbor distances, 
as shown in Figure 19. Note that such a histogram would 



Figure 18. Density probability of 75 spacecraft in family (5 R £ 
spacing) and single (or random) spacecraft distributions of orbit 
apogees, both averaged over random distributions in orbital phase, 
with all perigees at 3 R £ and apogees extending to 35 R £ . The 
probability peak at apogee results from slow motion there. 


— in the limit of a perfectly uniform grid of spacecraft (like 
that of a global simulation) — produce a delta function 
histogram with its peak density at the grid spacing dis- 
tance. We adopted the most probable nearest neighbor dis- 
tance as the effective spatial resolution of the constella- 
tion, and we sought a constellation whose distribution of 
nearest neighbor distances was as narrow as possible, and 
therefore as close as possible to an ideal grid of spacecraft. 
To avoid overly weighting the perigee legs of the orbits, 
nearest neighbor distances are computed only for the parts 
of the orbits with X<-10 R E . For 75 spacecraft spread over 
the range from 10 to 35 R E , in either distribution, the reso- 
lution is about 1 .2 R £ . Relaxing the requirement to 2.0 R E 
allows the number of spacecraft in the constellation to be 
reduced somewhat. When scaled to a distribution extend- 
ing from 7 to 40 R E , the mission science goal is a constella- 
tion of 100 spacecraft with resolution of 1.2 R £ . However 
the descope floor that just satisfies the mission require- 
ment for 2 R e mean sparing is 50 spacecraft. 

It can he seen from Figure 19 that as the number of space- 
craft decreases, the mean spacing increases, as expected. 
By distributing the apogees appropriately, it is possible to 
level the spacecraft probability, as suggested by Figure 
18, in the radial dimension. This produces a narrower dis- 
tribution of spacings, and a more uniform constellation dis- 
tribution overall, than shown here. 

3.1.2 Constellation Evolution 

One of the important considerations about this constella- 
tion configuration is how long it will stay together when 
the orbits evolve. To check this we have calculated orbits 
including moments of the Earths gravitational Field up to 
the fourth geopotential coefficient, J4, and including lu- 
nar and solar gravitational effects. 


26 


Magnetospheric Constellation Mission... 



Nearest Neighbor Distance (R E ) 


Figure 19. Histograms of nearest neighbor distances for a fleet of varying 
spacecraft with orbital apogees over the range from 7 40 R £ .. 


The evolution of the constellation orbits is affected by or- 
bited insertion uncertainties, by lunar and solar perturba- 
tions, and by differential precession among the different 
orbits. Insertion uncertainties will produce some initial 
spread of the orbits away from the targeted alignments. 
Lunar/solar perturbations will cause the outer orbits to 
drift away from the initial orbital plane. Figure 20 sum- 
marizes and illustrates the manner in which the orbit in- 
clination will evolve over the first 2 years. Differential pre- 
cession will spread the inner orbits relative to the outer 
orbits, splaying the overall constellation out in local time. 
Figure 21 illustrates and summarizes the orbital preces- 
sion effects for orbit inclination of 10° (and a slightly dif- 
ferent constellation design). 

3.1.3 Deployment Strategy 

The strategy for orbital deployment supports a single 
launch of up to 100 nanosatellites using a Dispenser Ship 
(DS) spacecraft that contains and deploys the individual 
nanosatellites. The DS is designed to fit within a Delta II 
fairing, and to hold up to 105 nanosatellites. The baseline 
mission design assumes that each nanosatellite is equipped 
with sufficient propulsion to reach its final orbit from a 
suitable drop-off orbit attained by the DS. The propulsion 
allocation between the DS and the nanosatellites is a sub- 
ject for future trades. 


3.1.4 Radiation Environment 

The expected radiation environment for spacecraft in the 
DRACO orbits is an additional consideration for the or- 
bital deployment strategy and has been evaluated. The total 
ionizing dose for a 2-year mission is 100 k-rad. Figure 22 
shows the radiation dose-depth curves for a 12.5 R £ apo- 
gee orbit. It can be seen that this environment, while sig- 
nificant, is comparable to other mission dose requirements, 
notably POLAR. This environment will require trades be- 
tween shielding (increased mass) and custom electronics 
(increased cost). 

3.2 Instrument Accommodations 

The ST-5 instrument provides the required mechanical ac- 
commodations for DRACO payload s footprints and fields 
of view. It also provides power and data storage. Table 4 
summarizes the instrument resource requirements. 

3.3 Spacecraft Description 

The information contained in this section is from NASA 
GSFC's ST-5 proposal, modified to accommodate the Mag 
netotail Constellation payload. The spacecraft platform is 
spin stabilized with body mounted solar arrays. The basic 
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Figure 20. Evolution of typical Magnetotail Constellation 
orbits over 2 years. 

structure is an octagonal prism 20.0 cm in height and 4 1 .6 
cm in diameter across the flat facets (44.0 cm including 
the stowed magnetometer) and 45.0 cm across the verti- 
ces. The overall height in the launch configuration with 
both the Global Positioning System (GPS) and whip-type 
Ultra High Frequency (UHF) antennas stowed is 29.0 cm 
from the base of the balance weights to the top of the X- 
band antenna. Figure 23 shows the deployed spacecraft. 
Table 5 summarizes the spacecraft component and mass 
list. Note that the ST-5 instrument mass is identical to the 


Magnetospheric Constellation instrument mass 
allocations. 

3.3.1 Subsystem Description and 
Performance 

Table 6 summarizes the predicted ST-5 subsystem sizing 
requirements, component capabilities, and integrated per- 
formance results for the spacecraft functions. Each sub- 
system is discussed in detail in the following. 

3. 3. 1.1 Spacecraft Mechanical Subsystem. The me- 
chanical subsystem consists of a structural bus, 
deployables, and a launch vehicle deployment/release 
mechanism. The graphite/cyanate ester composite bus 
structure consists of a removable, high-conductivity hon- 
eycomb sandwich top deck, to which most of the spacecraft 
equipment is mounted: 

• a composite honeycomb sandwich bottom deck; and 

• sidewalls constructed of thermally isolated filament 
wound composite shell. 

Components will be attached to the top and bottom decks 
with aerospace fasteners and/or adhesives. Aluminum ra- 
diation shielding, 0.218 cm, is provided for sensitive com- 
ponents Command and Data Handling (C&DH), transpon- 
der, power supply electronics, and A Formation Flying and 
Communications Instrument/Global Positioning System 
(AFFCI/GPS) to restrict the total dose to 50 k-rads. Soft- 
ware-commanded deployable structures include the mag- 
netometer boom and the GPS and UHF antennas. The mag 
netometer boom deployment is repeatable to 0.25 degrees 
(each axis). 



Figure 21. Drift of a constellation of 21 spacecraft over a 2-year period owing to differential precession of lines of apsides, 
for orbit inclination of UP. Spatial frequency of occurrence along each orbit is color coded in arbitrary units a) at launch , b) 
after 1 year, and c) after 2 years. The precession of the constellation shown here is much more rapid than for DRACO \ 
because this constellation has perigee at low Earth orbit. It is provided as an extreme case reflecting the fastest possible 
precession in the event that, lower perigees are used for DRACO. 
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Figure 22. Constellation DRACO radiation environment . Orbital 
parameters are given in terms of apogee and perigee altitudes , 
and correspond to 12.5 R £ x 3 R E geocentric distance. 


3. 3. 1.2 Electrical Power Subsystem. The Electrical 
Power Subsystem (EPS) consists of three major compo- 
nents: the solar array the battery, and the control elec- 
tronics. Triple junction Gallium Arsenide (GaA’s) solar cells 
with 26% beginning-of-life (BOL) efficiency and conduc- 
tive covers populate two separate cell strings on each of 
the eight solar panels (total array area of 832.0 cm 7 ). Each 
string is connected to a separate battery cell. Flexible cir- 
cuit interconnects provide electrical connection for the so- 
lar cells on one of the eight panels. Standard wire inter- 
connects are used on the remaining panels. 

A two-cell, 42 W-hr lithium ion battery supplies power dur- 
ing eclipse and for in-Sun "peaking" loads exceeding the 
solar array capability. Connecting each battery cell to a 
separate solar cell string incorporates charge protection. 

Rate solar cell string incorporates charge protection. String 
voltage is limited to 4.10V by solar array string regula- 


tors, which assure full charge at the maximum battery cell 
voltage of 4.10V and a maximum battery total voltage of 
8.2V. The C&DH, which monitors individual cell voltages 
and reduces spacecraft total load when necessary, accom- 
plishes discharge voltage control. Power conditioning and 
control electronics regulate solar array output and supply 
regulated power to spacecraft loads. 

3.3. 1.3 Thermal Control System. The thermal control 
system (TCS) is primarily passive and includes heaters, 
multilayer insulation (MLI) blankets, fillers, thermistors, 
and variable emittance coatings. The interior of the solar 
array panels is insulated with MLI as are the interiors of 
the top and bottom decks, aside from cutouts for a single 
radiator on each deck. The battery will be insulated to main- 
tain 25 °C at the minimum expected temperature of -30 
°C. One heater provides 0. 1 5W to the magnetometer and a 
second provides 0.5W to the transponder oscillator. Exter- 
nal coatings, including the radiator coatings, are electri- 
cally conductive, and all MLI layers are grounded to the 
structure. 

3.3. 1.4 Propulsion System. Trade studies were completed 
to explore the most efficient way to put each nanosatellite 
into its desired orbit. This study weighed the merits of 
putting all the propulsion on the deployer ship, all on the 
nanosatellites, or some combination of the two. Viable ver- 
sions of each scenario were judged to be nearly equal in 
merit. 

In case some propulsion is required for each nanosatellite, 
advanced engineering studies have already shown that 
such new systems are feasible. The nanosatellite propul- 
sion system consists of a miniature cold gas thruster sys- 
tem and a 15.2 cm diameter titanium Gaseous Nitrogen 
(GN2) propellant tank. Within the operating pressure range 
from 1500 psia BOL to 100 psia, the cold gas thruster sys- 
tem can provide a total velocity increment of 5.5 m/s. The 
system can deliver a minimum impulse bit of 45 mN-s at 
1000 psi. The system provides a maximum vacuum thrust 
of 1.32 N at 1500 psia and a maximum impulse of 109 N-s 
at a specific impulse of 60 seconds, providing a total delta 
Vee capability of 5.5 m/s. The ST-5 propulsion system will 
require modification if the individual satellites need to pro- 
vide the propulsion for apogee changes after deployment 
from the dispenser ship. 

3.3. 1.5 Guidance, Navigation, and Control Sub 
system. The Guidance, Navigation, and Control (GN&C) 
subsystem incorporates a mini +/- 90 degree field-of-view 
Sun sensor, a nutation damper, and also uses the science 
magnetometer. The attitude control system software actu- 
ates the cold gas thrusters via the C&DH subsystem to 
perform minor orbit adjustments and spin axis precession 
control. Spin-synchronized pulses are used for precession 
control, while continuous firing is used for orbit raising. A 
simple passive nutation damper damps thruster-induced 
nutation. 
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MISC MECH POWER C&DH THERMAL PROPULSION GN&C COMM INSTRUMENTS 


Mass Ea. Unit 


Quantity Component (kg) 


/ Magnetometer 0.150 

1 MAC Electronics, etc. 0.800 

l MAC Boom 0. 750 

7 Electron Detector 0.500 

7 Proton Detector 0.500 


1 Transponder, Xband 0. 1 50 

1 X- band Antenna 0.100 

1 Amp-Diplexer (HPA) 0.160 


1 Sun Sensor 0.250 

1 Nutation Damper 0.100 

1 CCNT 1.550 

1 GPS Antenna 0.100 

1 UHF Antenna 0.100 


Thruster 0. 130 

Tank 0.560 

Cold Gas 0.200 

Propulsion Misc 0. 5 00 


1 Thermal 0.760 


1 C&DH box 1-800 


/ Battery 1.000 

/ Solar Array 1.258 


1 Separation Device 0.200 

1 Balance Weights 0.500 

1 Structure 1.815 


7 Radiation Shielding 0.244 

7 Harnesses (all) 


Totals 


Total Mass 

Dimensions 


(kg) 

(cm) 

Mass of 

Subsystem Actual 

Subsystem 

0.150 

4x4x6 

2.700 

0.800 

10x20x5 


0.750 

50L x 3 


0.500 

10x10x7 


0.500 

10x10x7 


0.150 

5 x 7 x 5 x 5 

0.410 

0.100 

7L 


0.160 

2x4x4 


0.250 

2x2x2 

2.050 

0.100 

4x4x1 


1.550 

15x6x18 


0.100 

6x6x5 


0.100 

.2 x .2 x 20 


0.130 

0.560 

0.200 

0.500 

1.75 ODx 2.2 

1.440 

0.760 

N/A 

0.760 

1.800 

2 x 18 x 36 

1.800 

1.000 

1.258 

6.35 x 6.35 x 11.8 

2.258 

0.200 

0.500 

1.815 

N/A 

2.515 

2.684 

N/A 

4.013 

1.329 




17.946 


Table 5. ST 5 spacecraft component and mass summary. 
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Figure 23. Conceptual view of the Space Technology-5 
satellite for Magnetotail Constellation. 


3.3. 1.6 Command and Data Handling Subsystem. The 

C&DH subsystem is an Essential Servicing Node (ESN)- 
based system, which resides on one side of a double-sided 
card. The ESN multichip module (MCM) contains an ASIC 
chip as its central device. Embedded in the ASIC is a gen- 
eral purpose bus computer with the UTMC UT69R000 1 6- 
bit RISC processor as its core and a number of industry 
standards as its onboard peripherals. In addition to the 
devices contained in the ASIC, the MCM contains several 
memory storage devices used for instruction and datastor- 
age and an analog subcircuit used to support analog sig- 
nal measurements. 

The C&DH system sends commands, gathers telemetry 
data, stores data in Consulting Committee for Space Data 
Systems (CCSDS) packets, time tags data, and transmits 
stored data to the ground. The card consists of: 

• an ESN multichip module, acting as the MHz proces- 
sor which processes commands, makes data packets, 
and performs ACS computations; 

• multiple RS422 interfaces which gather telemetry and 
distribute commands to the science instruments and 
to the space-to-ground and space-to-space communica- 
tions systems; 

• a memory section containing 80 Mbits of D-RAM; 

• an ACS interface for the Sun sensor and thruster; 

• a housekeeping interface which reads analog data from 
the power system and thermistors and controls 
heaters; 


* a communications interface, which sends CCSDS packet 
telemetry to the ground and receives and processes 
CCSDS packet commands. 

3.3. 1.7 Communications Subsystem. The space-to- 
ground communications system consists of the Aero-Astro 
X-band transponder and diplexer. The transponder is used 
to downlink all science and housekeeping data to the 
ground once per orbit, to receive ground commands, and 
for two-way coherent Doppler tracking with the ground. 

3.3. 1.8 Flight Software Subsystem. The flight software 
is based on the operating system for the low-power ESN 
Central Processing Unit (CPU) used on the MAP mission. 
All flight software resides in a single 1 6-bit UT69R000 ESN 
microprocessor with a configuration similar to that of the 
MAP housekeeping remote services node, and uses the "C M 
language and assembly level programming. 

3.4 Dispenser Ship 

The details presented herein for the Magnetospheric Con- 
stellation DS are given by Lieberman [ 1 999b] . That report 
also references and draws upon the Magnetospheric Con- 
stellation Mission Document \Lieberman , 1999a]. The ac- 
tual mission design and DS designs have changed and may 
change further as a result of future studies. The orbits upon 
which the following details are based are not necessarily 
the orbits presented in Section 3.1. In case of discrepan- 
cies, the information given in earlier sections takes prece- 
dence over the materials in this section. 

3.4.1 Orbital Maneuvers of the Dispenser 
Ship 

Six methods for placing the nanosatellites in orbit were 
studied. The final orbit selected was the least complicated 
and minimized the nanosatellites propellant mass 
requirements. 

For the orbit selected, the Delta places the DS in an ellip- 
tical orbit whose parameters are perigee 185 km, apogee 
20 R £ with inclination of 28.7 deg. The DS performs an 
apogee burn and raises the perigee to 1000 km and makes 
a small change in the inclination. A second apogee burn 
raises the perigee to 3 R E and removes 21.2 deg of inclina- 
tion. The parking orbit of the DS is 3 R £ perigee, 20 R £ 
apogee, 7.5 deg inclination with the line of apsides lies in 
the ecliptic. From the parking orbit, several nanosatellites 
are released with their spin axes aligned to the DS perigee 
velocity vector. The nanosatellite kick motors are ignited 
at perigee. Some nanosatellites will raise their apogees and 
others will fire in retrograde direction to lower their 
perigees. 

Although the details presented in this report are based on 
delta-V capability for the nanosatellites, it has not been 
determined at this stage of formulation if this capability 
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1 Subsystem/Component 

Design 

Subsys tem/Compon en t 

Design 

Structure 


Propulsion 


Clearance Envelope (cm 3 ) 

43x31x45 

Velocity Increment (m/s) 

5.5 

Launch Loads (g’s each axis , flight) 
Min. Fundamental Axial Frequency 

>15 

Guidance. Navigation and Control 


>134 Hz 

• Pointing Accuracy, Sun Line to Spin Axis, Sun 

90+/ 3.0 

Min. Fundamental Lateral Frequency 
Radiation Total Dose Absorbed at 3 months 

134 Hz 
50 krads 

Acquisition Mode (deg.) 

• Spin Rate Knowledge, initial (% rpm) 

• Spin Axis Pointing Accuracy, Ground Based 

+/- 0.1 
2.5 



Control Mode +/- deg. 
* Knowledge (+/- deg.) 

LO 

Spacecraft Deployment 


Control (+/ deg) 

L5 

Separation Velocity (m/s) 

0.66 

* Spin Rate Knowledge On-Orbit (% rpm) 

+/ 0. 1 

Maximum Tlpoff Rate (deg. /sec) 

1.0 

• GPS Orbit Determination 3-D RMS (m) 

100 

Spin Rate at Separation (rpm) 

28+/- 

10% 

* PreDelta-V Spin Axis (Pointing +/- deg.) 

2.5 

Spin Rate Pre Delta V 

20+/- 

10% 





Command and Data Handling 


Power (Solar Array at 25° C) 


* PROM (kWord) 

32 

Power required „ orbit average (W) 

21.4 

* EEPROM (kWord) 

128 

Energy Storage Required (W hr) 

42.0 

* Data RAM (kWord) 

64 

Voltage Regulation Range 

5.0 +/- 

• Instruction RAM (kWord) 

164 


4.0% 

• Data Recorder Size (Mbits) 

80 

Thermal (Solar Flux - 1285 w/m 2 Cold \ 




1418 w/m z Hnl 




Internal Equipment Operating (C) 

5 to +25 

Communications (Space Ground) 

100 

Internal Equipment Non-Op (C) 

5 to +28 
77 to +28 

Data Rate (at 5.0 Re, 1.5 W, BER lx 10- 5 (kbits/ 
sec) 

5.5 

Solar Array Operating °C) 

Link Margin with 1 1 M Antenna (dB) 

Solar Array Non Op (C) 

77 to +40 

Total Data Downlinked (Mbit) 

640 

Duration of Elipse (hrs.) 

2.1 

Space to Ground Transmit Time (min.) 

107 

Radiator Area for 25 W (cm2) 

600 




Table 6. Spacecraft requirements and design capabilities 


will be the final nanosatellites requirement. Another op- 
tion under consideration is to design the DS with the ca- 
pability to ensure proper nanosatellite orbit insertion. 
Having the nanosatellite orbits determined solely by how 
and when they are released from the mother ship could 
have significant impact on the requirements and design 
for a release mechanism. 

3.4.2 Nanosatellite Orbital Maneuvers and 
Mission Operations 

Four nanosatellite release mechanisms have been proposed. 
The kinematics of release were simulated; the simulations 
show that all four will provide a noninterference separa- 
tion of the nanosatellite from the DS. All four require a 


discrete electrical command, which will cause the force 
holding the nanosatellite to the DS to cease to exist in a 
short period of time — approximately 1 msec. The com- 
mand will be given simultaneously to up to 14 
nanosatellites. When released, the rotating nanosatellites 
will "fly” tangentially away from the DS. The angular rate 
will be approximately the same as that of the DS. The com- 
mand receivers of the separated units will be turned on; 
however, only one will have its transmitter on. After the 
nanosatellite's kick motor has fired (the one with the ac- 
tive transmitter), the ground station — using the nominal 
burn trajectory for a reference orbit — will point its an- 
tenna at the nanosatellite. The Doppler shift of the trans- 
mitted signal will be used to determine the nanosatellite 
orbital velocity and by analysis the ground station will be 
able to determine the nanosatellites orbit. 
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When the nanosatellites are dispensed, they will be se- 
lected such that no two satellites will be transferred into 
an orbit with same apogee. This requirement will assure 
1) the ground station has adequate time to track the 
nanosatellites and, 2) the true anomaly of satellites in the 
same orbit will be different. If 10 nanosatellites are dis- 
pensed at 2 R e separation, no more than four will be in an 
orbit whose apogee is less than 20 R £ . Thus, the ground 
station will be able to get the 30 minutes of tracking per 
nanosatellite needed to get “good” data for each of the four 
“lower” apogee nanosatellites before they leave the search 
field. The ground station can then turn its attention to the 
remaining six nanosatellites destined for the higher apo- 
gees. Approximately 200 (325-120) minutes (33 minutes 
per satellite) are available to track the remaining six 
nanosatellites. 

3,4.3 Dispenser Ship Overview 

The heart of the system is the C&DH system. The X-band 
receiver, located in the Communications Subsystem, re- 
ceives the ground commands. They are decoded and stored 
in the C&DH as relative or absolute time commands. The 
ACS sensors send attitude signals to the C&DH, where 
they are processed to determine the actuator commands 
(five thrusters) that are needed to maintain attitude con- 
trol or to maneuver the DS to a new attitude. The C&DH 
issues release commands to the nanosatellites via the arm 
and fire relays. The C&DH also sends commands to the 
main thruster during each of the two orbit change maneu- 
vers. Finally, the C&DH processes and formats housekeep- 
ing, status, and health and safety telemetry for transmis- 
sion to the ground. Each subsystem is briefly described in 
this section; a more detailed discussion of the dispenser 
ship can be found in the full report \Lieberman, 1999a|. 

3.4.3. 1 Mass and Power Summary. Table 7 summa- 
rizes the mass and power budgets of the DS concept. The 
mass and power figures apply to the DS exclusive of the 
DRACO nanosatellites. 

3.4.3.2 Mechanical Configuration. Figure 24 shows the 
mechanical configuration of the nanosatellite DS.The blue 
ring at the top is the 130 watt solar array. Above the array 
is the X-band omni antenna. Located on the top ring are 
three coarse Sun sensors, Earth sensors and a fine Sun 
sensor. 

In this version, 92 nanosatellites are located in individual 
bays. There are 14 bays per tier and seven tiers, a total of 
98 bays. Since the lift mass of the Delta only permits a 
payload of 92 nanosatellites, six bays will be empty. Lo- 
cated in the area around the base of the DS are the main 
thruster, three additional coarse Sun sensors, and four 
small thrusters. The propulsion tanks, power regulators 
and the avionics are housed in the interior of the DS. 


3. 4.3.3 Attitude Control Subsystem (ACS) 
Description. The ACS uses an Earth sensor and fine Sun 
sensor for attitude references and to measure the spin rate. 
The ring laser gyro is used as attitude reference to precess 
the spin axis to a desired attitude and for sensing and con- 
trolling the nutation of the DS. The normal mission atti- 
tude of the spin axis is perpendicular to the ecliptic except 
during the portion of the orbit when the nanosatellites are 
to be released. Before release, the DS will be commanded 
to precess the spin axis to be parallel to the velocity vector 
at perigee. The nanosatellites will be released close to peri- 
gee. Precession will be accomplished by actuating the bi- 
propellant thrusters. 

The six coarse Sun sensors’ output provides the position of 
the Sun with respect to the DS. This information from the 
coarse Sun sensors will be used during the initial stabili- 
zation phase after the DS has been released from the launch 
vehicle. To precess the spin axis of the DS from the release 
attitude to the normal operational attitude (perpendicu- 
lar to the ecliptic), the thrusters are used. 

3. 4. 3.4 Command and Data Handling (C&DH) 
Subsystem. The C&DH avionics module is based on a cen- 
tral processor. The subsystem also has a solid-state recorder 
to store data that will be transmitted to the ground sta- 
tions for distribution to the end user. The avionics module 
includes interfaces to the ACS, the power subsystem, and 
the communications subsystems as well as the release ac- 
tuators. The avionics module sends commands, gathers te- 
lemetry data, processes the data, stores the data using 
CCSDS packets, time tags the data, and transmits the 
stored data to the ground. The avionics module is a proces- 
sor-based system which processes commands, packetizes 
the data, performs the timekeeping functions, performs 
ACS computations and controls power The avionics mod- 
ule contains 2.0 Gbits of Random Access Memory (RAM), 
and interfaces to the ACS subsystems Earth/Sun Sensor 
(analog data), Attitude Thruster (pulse command). Orbit 
Thruster (single pulse), the latch valves and the pyro ac- 
tuated valves. The housekeeping functions include read- 
ing the thermistors, pressure transducers, and voltage and 
current monitors, and gathering other analog telemetry. 
These functions will be handled using a 12-bit A/D con- 
verter. The Uplink/Downlink interface is connected directly 
to the transmitter and receiver. This interface sends CCSDS 
packet telemetry to the ground antennae, and receives 
CCSDS packet commands from the ground antennae and 
processes them. 

The avionics module is an MCM stack. It is composed of 
six functional interfaces, each of which performs separate 
functions required by the avionics module. These functions 
are the processor, the memory, the actuator interface, the 
ACS interface, the housekeeping interface, and the com- 
munications interface. The development of these interfaces 
will require the extensive use of Field Programmable Gate 
Arrays (FPG As). These FPGAs will be developed in a low- 
power technology. 
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Subsystem 

Mass 

(kg) 

Power 

(watts) 

ACS 

7.1 

13.8 

C&DH 

12 

2 

Comm unications 

7.8 

7.2 

Power 

5.6 

32 

Propulsion 

50 

- 

Harness 

4.5 

0 

Thermal 

4.5 

5 

Structure 

131 

0 

Total as Designed 

216.9 

60 


Table 7. DS final as- designed mass and power. 


3. 4. 3. 5 Communications Subsystem. The communication 
system for the DS consists of a transmitter, receiver and an 
omni antenna. It operates in the X band. The link analysis 
shows positive margin with an 1 1 -meter ground station an- 
tenna and an omni antenna on the DS. The C/T ratio for the 
ground segment must be 35.4 dB/K in the X band. The Effec- 
tive Isotrophic Radiated Power (EIRP) from the DS is 1.76 
dbw (1.5 watts to the antenna). To be compliant with the 
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Figure 24. Conceptual mechanical configuration for nanosatellite 
dispenser ship. 


system requirement, using one-way Doppler for track- 
ing, the oscillator must be stable to 5 x 10* per day. The 
down link transmission rate can be as high as 100 Kb/ 
sec and the command rate is baselined as 1 Kb/sec. (NB: 
downlink rates of up to 640 Kbps have been determined 
to be feasible, since this study). 

3.4. 3.6 Power Subsystem. The power subsystem con- 
sists of three main elements: the solar array which con- 
verts solar energy to electrical power, the regulators, 
which, as the name implies, reduce the variation of the 
voltage on the power bus, and the lithium hydride bat- 
tery. The output to the bus will be 28 +/- 7V DC. The 
solar array is sized to deliver 130 watts off the array at 
BOL. The power subsystem provides no power to the 
nanosatellites. The power subsystem provides power 
through eclipses up to 40 minutes. 

3.4.3. 7 Propulsion Subsystem. The DS propulsion pro- 
vides 202 kg of bipropellant and operates while the DS 
is spinning at 30 rpm. The system will be used for add- 
ing delta- V, preccssing the spin axis, controlling the spin 
rate and maintaining attitude control. The large thruster 
will be used for delta-V maneuvers and the smaller 
thrusters will be used for the other functions. The pyro 
valves are opened after the DS has been placed in orbit 
by the launch vehicle. The thrusters will not fire until 
the latch valves are commanded open. The main engine 
latch valve will be kept closed and opened when a delta- 
V maneuver is imminent. The other latch valves will be 
left open and closed only if a thruster malfunctions. The 
thrusters are fired infrequently and therefore the orbital 
average power is negligible. 

3. 4. 3. 8 Structure. The DS structure consists of the in 
ner structure and seven berthing rings. The inverted cone 
at the base of the internal support transition cylinder 
interfaces with the launch vehicle payload adapter fit- 
ting. Attached to this cylinder are seven vertical gus- 
sets, which support the seven nanosatellite berthing 
rings. Detailed views of these structures are shown in 
the full report [Lieberman, 1999a]. 

3. 4. 3. 9 Inner Structure. The propulsion system’s four 
tanks and the large thruster are carried by the inner 
structure. The four small thrusters are located on the 
lower berthing ring. Two of the thrusters are for spin 
control, a third is used for precession control and the 
fourth is a backup for the precession control thruster. In 
addition, the inner structure carries the three coarse Sun 
sensors. 

3.4.3.10 Thermal Subsystem. The DS thermal require- 
ments can be accommodated with a passive thermal con- 
trol system (blankets, heaters, and coatings). The MLI 
will be used on top and bottom to reduce losses and gains 
from the environment. The cylindrical sides (shell) of the 
DS are not insulated. Heaters and MLI may be used to 
maintain propulsion lines and thrusters during orbit 
changes. 
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A top-level thermal analysis has been completed. For com- 
ponents that are mounted externally, have a high watt den- 
sity or have tight temperature control requirements, a more 
detailed model would need to be completed. 

Preliminary results indicate that the nanosatellites 
mounted in the top ring run slightly warmer than the oth- 
ers, while the ones in the bottom ring run the coolest. The 
DS temperature is set by the high energy absorbed by the 
shell and its conductive coupling to the birdcage. The de- 
sign described is reliable and appears to meet tempera- 
ture requirements for the DS as well as the nanosatellites. 
Further parametric analysis will be conducted with a more 
detailed model 

3.5 Integration and Test 

Of particular importance to the DRACO mission will be 
innovation in the way in which the payloads 
(nanosatellites) are integrated to the dispenser ship. The 
normal way of integrating payloads to the main spacecraft 
would result in an unacceptably long and costly integra- 
tion period. This section discusses w 7 ays of approaching this 
issue. It will be assumed that each nanosatellite bus and 
instrument suite will be delivered as a single, ready-to- 
integrate system. Furthermore, each nanosatellite must 
not take longer than 1 or 2 days to be fully integrated to 
the DS. This is necessary because the integration of 92 
nanosatellites taking 2 days each would be 184 days, or 6 
months of 7-day per week operation if the nanosatellites 
are integrated in a serial fashion. If several nanosatellites 
could be integrated in parallel, this would reduce the inte- 
gration time, but could increase integration cost due to 
higher manpower requirements. 

One approach to mitigating this issue is to design the 
nanosatellite for rapid integration thru the use of highly 
automated processes and test equipment. This would re- 
quire development of automated test equipment and pro- 
cesses, or would require that the design process take into 
account the use of existing test equipment and processes 
to optimize the ability to quickly test the "sciencecraft.” 

Another approach is to accept more risk in the process by 
only completely testing every n th unit and performing re- 
duced testing on the other units. This approach means that 
more risk must be accepted. 

3.6 Mission Operations 

The operation of a constellation of identical nanosatellites 
requires different concepts than for single spacecraft mis- 
sions. The nanosatellites may be constrained in the amount 
of functions they can perform on board by the lack of pro- 
cessing power. The large number of spacecraft requires au- 
tomation on the ground in order to keep the staffing to a 
reasonable level and offers an opportunity for changes in 
risk management. The loss of a few spacecraft over the 
mission lifetime is tolerable for this mission. 


Routine operations commence after the nanosatellites have 
reached their operational orbit. Anomalies or infrequent 
events such as eclipse only interrupt routine operations. 

3.6.1 Space/Ground Communications 

There are several options for space/ground communications 
services: dedicated stations, commercial networks, NASA 
stations, large antennas, smaller antennas, etc. This op- 
erations concept assumes the use of a commercial network 
of 1 1 -meter antennas. The commercial network would have 
a number of these antennas distributed around the world. 
Communications would only occur around perigee. At times 
during the mission, there may be as many as two dozen 
spacecraft near perigee at the same time. The multiple sta- 
tions would be able to handle some of these spacecraft. 
Spacecraft with smaller apogees could skip a contact, and 
dump their data at the next perigee. Loss of some data 
during these infrequent "traffic jams’' would be acceptable 
— the overall data completeness requirement will be - 95%. 

The ground stations will provide telemetry commanding, 
and tracking services. The real-time housekeeping data will 
be extracted from the downlink data and sent to the op- 
erations center in real-time. Playback science and house- 
keeping data and tracking data will be sent to the opera- 
tions center after the contact, sharing the bandwidth with 
other station users. 

3.6.2 Planning and Scheduling 

The scheduling of the ground stations depends on how the 
ground station services are provided. One option would be 
to let the ground station provider schedule the contacts 
within provided guidelines for data completeness and fre- 
quency of tracking data collection. The ground station pro- 
vider would have the flexibility to adjust the contacts within 
these guidelines, which potentially could result in lower 
costs to the project. The ground station provider would 
provide the schedule to the operations team so that they 
could identify uplink opportunities and support data 
accounting. 

3.6.3 Commanding 

The spacecraft will be commanded every contact to ini- 
tiate the downlink of data. Since the ground station sched- 
ule is subject to change, particularly for the spacecraft with 
distant apogees, the spacecraft cannot initiate the down- 
link. The instrument scripts may be updated once every 
month. The instrument script may be as simple as timed 
tagged commands or could he a more sophisticated set of 
instructions for the onboard system to interpret. Software 
loads may be performed if problems are identified after 
launch that can be addressed by the flight software. The 
data system will automatically upli nk new loads or instru- 
ment scripts to the constellation. The loads will have a time 
window within which each spacecraft must be updated. 
The system will automatically send and verify the load 
during regularly scheduled contacts. The operators will be 
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alerted of any spacecraft that is not updated when the win- 
dow closes. 

3.6.4 Housekeeping Data Processing 

The data system will automatically process the housekeep- 
ing data and identify limits violations and improper con- 
figurations. The system will be capable of processing data 
from multiple spacecraft simultaneously. The data system 
will also process status information from the ground sta- 
tion and identify any missing or poor quality data. Prob- 
lems will be ranked according to severity. Significant prob- 
lems will result in an immediate alert to the operations 
team. If the operations team members are not present, the 
system will use a pager to notify them. The data system 
will process all of the housekeeping data and produce stan- 
dard analysis reports. These reports will compare telem- 
etry parameters across time and across spacecraft. Space- 
craft engineers will be able to generate custom reports as 
well. 

3.6.5 Science Data Processing 

The science data is automatically level -zero processed when 
it is received from the ground station. This processing re- 
moves overlaps in the data and identified missing data. 
Operators are alerted in the event that a significant amount 
of data is missing. The system will have data accounting 
tools to assist the operators in visualizing the data loss 
and its distribution over the constellation. The operators 
will work with the ground station network to attempt to 
recover lost data that exceeds requirements. They may 
adjust the scheduling priorities if some part of the constel- 
lation is losing more data than other parts. The data will 
be automatically processed into standard products and 
archived in the science center. It will be distributed to sci- 
ence users electronically or, for large amounts of data, on 
physical media. 

3.6.6 Operations Staffing 

The mission operations staff (not including the science op- 
erations staff) for the Magnetospheric Constellation 
DRACO can be relatively small. The spacecraft are per- 
forming survey missions and have limited reconfiguration 
capability. The ground systems are automated, perform- 
ing all routine functions. The operations team will prima- 
rily handle exceptions that have been flagged by the data 
system. The operations team’s primary concern will be for 
common problems that have the potential to affect all of 
the members of the constellation. The nominal operations 
staff will be about six people. It will include one or two 
spacecraft engineers, a science liaison member, a data ac- 
counting person, and an orbit expert. The operations staff 
will work 40 hours per week, and at least one member will 
be on call for automated alerts during off-hours. 


3.6.7 Special Operations 

At launch, the operations team will be augmented with 
personnel from the spacecraft development and test func- 
tions. The ground data system may include more ground 
stations or larger ground stations, to check out and track 
the spacecraft after deployment. The ground system will 
also have to operate the DS during this phase of the mis- 
sion. The deployment operations may be adjusted (e.g., to 
change the number of spacecraft deployed at a time, or 
change the frequency of deployment), based on the experi- 
ence with the first few deployments. 

In the event of an anomaly, the operations team may re- 
quest the assistance of the spacecraft developers to iden- 
tify the cause of a problem and to develop corrective action 
or workarounds. The spacecraft developer will be respon- 
sible for maintaining the flight software throughout the 
mission. To assist in the anomaly resolution, the opera- 
tions team may use larger ground antennas to communi- 
cate with the spacecraft over longer periods of time, rather 
than just at perigee. 

3.7 Status and Schedule 

The current STP funding profile supports a launch in Sep- 
tember 2010. Implementation will begin in 2006. Technol- 
ogy development Is expected to continue through 2005. The 
schedule for the project is shown in Figure 25. 

4.0 Mission-Enabling Technology 
Development 

To realize the Magnetospheric Constellation mission, a 
number of enabling technologies must be developed and 
matured. These technologies are listed below: 

• Low-mass, low-power, C&DH systems; 

• Miniaturized, low-power, radiation-hardened analog 
and digital electronics; 

• Advanced propulsion systems; 

• Avionics on a chip; 

• Lightweight, high-capacity batteries; 

■ Lightweight, low-power communications components; 

• Software tools for autonomous ground operations of 
satellite constellations; 

• Data acquisition, distribution and visualization from 
constellations. 

In addition to these technologies, new tools and method- 
ologies must be developed for manufacturing large num- 
bers of research grade instruments and satellites at low 
cost, and testing and integrating these large numbers of 
nanosatellites in a reasonable amount of time. These par- 
ticular issues will be discussed in greater detail in the sec- 
tions to follow. 
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4.1 Instrumentation 

Development of instrument technologies and instruments 
of mass and power commensurate with a nanosatellite class 
mission are well underway by Department of Defense 
(DOD). NASA and various other institutions and agencies, 
both U.S. and international. Some instruments, such as 
magnetometers, are already close to the mass and power 
needed for a constellation mission. The issues of most sig- 
nificance in the development of instruments for a constel- 
lation mission are recurring cost and capability to be manu- 
factured in large quantities in a reasonable period of time. 
The current paradigm of a 2-3 year instrument develop- 
ment, fabrication and qualification period with costs typi- 
cally in the several hundred thousand to over $1 million 
cost range is clearly unacceptable when one wishes to fab- 
ricate and qualify as many as 100 identical sets of instru- 
ments. It is this "manufacturability” issue that most needs 
to be addressed by the Magnetospheric Constellation 
DRACO mission, especially in light of the possibility that 
some instrument technologies may require new designs to 
achieve the cost and "manufacturability" requirements. 

The strategy to be pursued by the Magnetospheric Con- 
stellation mission includes the following options: 

• NASA Research Announcements (NRAs) initiated and 
funded through DRACO formulation efforts; 

• Leveraging of NASA Headquarters-initiated and 
funded NRAs (directed DRACO funding to supplement 
these efforts); 

• Participation in, or leveraging of. other proposal oppor- 
tunities such as the Cross-Enterprise Program; 

• Spinoff from the ST-5 nanosatellite Constellation T rail- 
blazer project under the New Millennium Program; 

• Restructuring the typical spending plan for STP mis- 
sions in order to fund preliminary designs for early 
mission development. 

Another approach under consideration involves the 
partnering of scientists and instrument developers with 
the developers and builders of spacecraft buses to create 
an integrated "sciencecraft.” This approach will be dis- 
cussed in a little more detail in the next section. 

4.2 Nanosatellite Development 

In contrast to instruments, a number of technologies (as 
listed above) needs to be developed and matured before a 
nanosatellite meeting the DRACO requirements can be re- 
alized. As with the instrument technologies, a number of 
agencies and institutions, both domestic and international, 
are actively working on nanosatellite technologies. In many 
cases, these technologies are being funded at the subsystem 
level of development (e.g., avionics on a chip or miniatur- 
ized propulsion systems), with emphasis placed on the tech- 
nology of most interest by the organization providing the 
money. The DRACO mission will track these activities and 
wherever possible leverage these development efforts. One 


potential near-term opportunity is to work closely with the 
recently awarded ST-5 mission, nanosatellite constellation 
Trailblazer. Efforts are cilready underway to establish a 
relationship with the ST-5 project office at Goddard Space 
Flight Center. 

The same issues of nonrecurring cost and “manufactur- 
ability” apply to the spacecraft. In addition, however, there 
is the issue of integration and testing at the laboratory 
level (instruments and spacecraft). Although some efforts 
at development have begun in these areas (e.g., National 
Research Opportunity (NRO) Directors Innovation Initia- 
tive NRO000-99-R-0176), the issue as it applies to DRACO 
must be addressed. Although the spacecraft "manufactur- 
ability” issues can be addressed through the use of NRAs 
and Requests for Proposals (RFP’s), a more integrated ap- 
proach has been suggested for the DRACO mission. 

It is thought that all of the manufacturability issues might 
be addressed by teaming a science group with a spacecraft 
vendor. In this fashion, an integrated approach can be taken 
from the beginning, at the design level. Innovative meth- 
ods of testing must be developed to greatly reduce the time 
required; this can be accomplished through hardware de- 
sign, software development or a combination of both. Test 
methods and hardware can be designed concurrently with 
the instrument and spacecraft development process. In- 
strument and spacecraft designs would be influenced by 
how easily and quickly integration could be done. At the 
present time, it is not clear what mechanism would be used 
to create or encourage the partnership of the science teams 
and the spacecraft teams. 

4.3 Constellation Operations Management 
and Autonomy 

In order to manage a large constellation of nanosatellites 
during the operations phase at a reasonable cost, a signifi- 
cant degree of autonomy on the ground as well as on board 
the nanosatellites will be required. There is already an ex- 
isting base of experience and knowledge, along with the 
hardware and software tools, that resides in the commer- 
cial world — in particular, the Iridium communications 
constellation. It is expected that the DRACO formulation 
efforts will learn and build upon this knowledge base, de- 
veloping new tools and methods as required for a scien- 
tific constellation. To the extent that the commercial world 
is aware of the push for future science missions that con- 
sist of large constellations of spacecraft, there is already 
some evolution of the ground operations tools and meth- 
odologies to accommodate this. 

In the area of spacecraft autonomy, there is already a sig- 
nificant amount of interest in this area, with work being 
funded through numerous channels. The DRACO mission 
will require as a minimum the capability for onboard er- 
ror detection and correction and the ability to operate for 
days or weeks at a time without ground intervention, but 
to a large extent these capabilities already exist. A greater 
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Figure 25. Magnetotail Constellation development schedule. 


degree of sophistication will be required for the DRACO 
mission and it is possible that new requirements will be 
levied in this area as the formulation effort progresses. At 
the current time, however, the exact degree of autonomy 
and detailed requirements for the DRACO Mission has not 
been determined beyond the basic requirements mentioned 
above. 

4.4 Advanced Data Assimilation and 
Handling Techniques 

It is widely recognized that the data sets generated by con- 
stellation missions require data management and handling 
techniques that are outside the current experience with 
space missions. Nevertheless, the problem is not qualita- 
tively different from that associated with the conduct of 
atmospheric meteorology research and monitoring. 

The most fundamental problem is the visualization of both 
scalar and vector field information from a large number of 
observing stations. Some techniques developed for meteo- 
rology will be directly transferable to use with constella- 
tion data sets. Others will require new techniques, for ex- 
ample, the visualization of magnetic field information that 
is not present in atmospheric meteorology. 


The next step beyond basic visualization is the assimila- 
tion of the data into global circulation models. This will 
allow for instantaneous comparison with such models, ex- 
trapolation of the observations into other parts of the sys- 
tem outside the constellation, and projection or forecast- 
ing of system behavior for comparison with actual observed 
system evolution. This is a complex process, which has been 
discussed in much more detail in the section on science 
objectives, above. 

Because the treatment of constellation data is so foreign 
to space physics, there is a significant need for the devel- 
opment of appropriate tools, borrowing wherever possible 
from prior work along the same lines that has already been 
completed in meteorology. As an integral part of the tech- 
nology development effort for DRACO, developments along 
these lines should be supported at appropriate institutions 
with expertise in visualization and assimilation of data. 
Researchers who are active in the global simulation of 
magnetotail processes should undertake a significant 
amount of this work. They are in a position to generate 
simulated data sets, to introduce realistic errors in them, 
and then to develop the means to reassimilate the simu- 
lated data back into the same model that generated it, or 
other independent models. 


38 


Magnetospheric Constellation Mission... 






































Appendix A: Selected References 

Aggson. T.L., J.P. Heppner. and N.C. Maynard, Observa- 
tions oflarge magnetospheric electric fields during the 
onset phase of a substorm, J. Geophys. Res., 88, p. 3981 , 
1983. 

Akasofu, S.-I.,The development of the auroral substorm, 
Planet. Space Sci., 12. pp. 273-282, 1964. 

Angelopoulos, V., C.F. Kennel, F.V. Coroniti, R. Pellat, H E. 
Spence, M.G. Kivelson, R.J. Walker, W. Baumjohann, 
W.C. Feldman, J.T. Gosling, and C.T. Russell, Charac- 
teristics of ion flow in the quiet state of the inner plasma 
sheet, Geophys. Res. Lett., 20, pp. 171 1-1714, 1993. 
Angelopoulos, V., C.F Kennel, F.V. Coroniti, R. Pellat, M.G. 
Kivelson, R.J. Walker, C.T. Russell, W. Baumjohann, W.C. 
Feldman, and J.T. Gosling, Statistical characteristics 
of bursty bulk flow events, J. Geophys. Res., 99, pp. 
21257-21280, 1994. 

Angelopoulos, V., F.V. Coroniti, C.F. Kennel, M.G. Kivelson, 
R.J. Walker, C.T. Russell, R.L. McPherron, E. Sanchez, 
C.T. Meng, W. Baumjohann, G.D. Reeves, R.D. Belian, 
N. Sato, E. Friis-Christensen, PR. Sutcliffe, K. Yumoto, 
and T. Harris, Multi-point analysis of a bursty bulk flow 
event on April 1 1, 1985, J. Geophys. Res., 101 , pp. 4967- 
4989, 1996. 

Angelopoulos, V., T.D. Phan, D.E. Larson, F.S. Mozer, R.P. 
Lin, K. Tsuruda, H. Hayakawa, T. Mukai, S. Kokubun, 
T. Yamamoto, D.J. Williams, R.W. McEntire, R.P. 
Lepping, G.K. Parks, M. Brittnacher, G. Germany, J. 
Spann, H.J. Singer, and K. Yumoto, Magnetotail flow' 
bursts: Association to global magnetospheric circula- 
tion, relationship to ionospheric activity and direct evi- 
dence for localization, Geophys. Res. Lett., 18, p. 227 1 , 
1997. 

Ashour-Abdalla, M., L.M. Zelenyi, V. Peroomian, and R.L. 
Richard, Consequences of magnetotail ion dynamics, J. 
Geophys. Res., 99. pp. 14891-14916, 1994. 
Ashour-Abdalla, M., L.M. Zelenyi, V. Peroomian, R.L. Ri- 
chard, and J.M. Bosqued, The mosaic structure of 
plasma bulk flows in the Earth’s magnetotail, J. 
Geophys. Res., 100, p. 19191, 1995. 

Baker, D.N.,T.I. Pulkkinen, R.L. McPherron. J.D. Craven, 
L.A. Frank, R.D. Elphinstone, J.S. Murphree, J.F. 
Fennell, R.E. Lopez, and T. Nagai, CD AW 9 analysis of 
magnetospheric events on May 3, 1986: Event C, J. 
Geophys. Res., 98, pp. 3815-3834, 1993. 

Baker, D.N., T.I. Pulkkinen, V. Angelopoulos, W. 
Baumjohann, and R.L. McPherron, Neutral line model 
of substorms: Past results and present view 7 , J. Geophys. 
Res., 101. p. 12975, 1996 

Baumjohann, W., and G. Paschmann, Determination of the 
polytropic index in the plasma sheet, Geophys. Res. 
Lett., 16, pp. 295-298, 1989. 

Baumjohann, W. p G.C. Paschmann, and C.A. Cattell, Aver- 
age plasma properties in the central plasma sheet, J. 
Geophys. Res., 94. pp. 6597-6606, 1989. 

Baumjohann, W., G.G. Paschmann, and H. Luhr, Charac- 
teristics of high-speed ion flows in the plasma sheet, J. 
Geophys. Res., 95, pp. 3801-3809, 1990. 


Baumjohann, W., A short update on the present-day near- 
Earth neutral line model, J. Atmos. Sol. -Terr. Phys., 62, 
2000 . 

Berchem, J., J. Raeder, and M. Ashour-Abdalla, Interac- 
tive visualization of numerical simulation results: A tool 
for mission planning and data analysis, in: Visualiza- 
tion techniques in space and atmospheric sciences, E.P. 
Szuszczewics and J.H. Bredekamp Eds. NASA SP-519, 
1995. 

Birn, J., M. Hesse, and K. Schindler, MHD simulations of 
magnetotail dynamics, J. Geophys. Res., 101 (A6), 
p. 12939, 1996. 

Borovsky, J.E., R.C. Elphic, H.O. Funsten, and M.F. 
Thomsen, The Earths plasma sheet as a laboratory for 
flow turbulence in high-beta MHD, J. Plasma Phys., 
57, p. 1, 1997. 

Calder, Nigel, The Weather Machine, Penguin Books, 
Middlesex, UK, 1974. 

Cattell, C., I. Roth, and M. Linton, The effects of low fre- 
quency waves on ion trajectories in the Earth’s 
magnetotail, Geophys. Res. Lett., 22, pp. 3445-3448, 
1995. 

Dungey, J., Inaugural Lecture as Professor of Physics at 
Imperial College, 1966: quote courtesy ofW. J. Hughes. 

Fairfield et al., Geotail observations of substorm onset in 
the inner magnetotail, J. Geophys. Res., 103, p. 103, 
1998. 

Frank, L.A., and W.R. Paterson, Survey of electron and ion 
bulk flows in the distant magnetotail with the Geotail 
spacecraft, Geophys. Res. Lett., 21, pp. 2963-2966, 1 994. 

Ghil, M. and P. Malanotte-Rizzoli, Data assimilation in me- 
teorology and oceanography, Adv. Geophys., 33, p. 141, 
1997. 

Ghil, M., K. Ide, A.F. Bennett, P. Courtier, M. Kimoto and 
N. Sato (Eds.), Data Assimilation in Meteorology and 
Oceanography: Theory and Practice, Meteorological So- 
ciety of Japan and Universal Academy Press, Tokyo, 
496 pp.,1997. 

Goertz, C.K., and W. Baumjohann, On the thermodynam- 
ics of the plasma sheet, J. Geophys. Res., 96, pp. 20991 - 
20998,1991. 

Goodrich, C.C., et al., Interaction between the inner and 
mid-tail regions in MHD simulations of substorms, EOS 
Trans. AGU, 80(46), p. 871, 1999. 

Haerendel, G., Disruption, ballooning or auroral avalanche 
on the cause of substorms, in First International Con- 
ference on Substorms, ESA SP-335, p. 417, 1992. 

Hayakawa, H., A. Nishida, E.W. Hones, and S.J. Bame, Sta- 
tistical characteristics of plasma flow in the magneto- 
tail, J. Geophys. Res., 87, pp. 277-283, 1982. 

Henderson, M.G., G.D. Reeves, H E. Spence, R.B. Sheldon, 
A.M. Jorgensen, J.B. Blake, and J.F. Fennell, First en- 
ergetic neutral atom images from POLAR, Geophys. 
Res. Lett., 24, p. 1167, 1997. 

Hones, E.W., Jr., The magnetotail: Its generation and dis- 
sipation, Physics of Solar Planetary Environments, ed. 
by D.J. Williams, 557-57 1, AGU, Washington DC, 1976. 

Hones, E.W., Jr., G. Paschmann, S.J. Bame, J.R. Asbridge, N. 
Sckopke, and K. Schindler, Vortices in magnetospheric 
plasma flow, Geophys. Res. Lett., 5. p. 2069, 1978. 


Magnetospheric Constellation Mission... 


39 


Hones, E.W., Jr., and K. Schindler, Magnetotail plasma 
flow during substorms: A survey with IMP 6 and IMP 
8, J. Geophys. Res., 84, pp. 7155-7169, 1979. 

Hones, E.W., Jr., J. Birn, S.J. Bame, J.R. Asbridge. G. 
Paschmann, N. Sckopke, and G. Haerendel, Further 
determination of the characteristics of magneto- 
spheric plasma vortices with ISEE 1 and 2, J. 
Geophys. Res., 86 . p. 814, 1981 . 

Hoshino, M.,T. Mukai, A. Nishida, Y. Saito, T. Yamamoto 
and S. Kokubun, Evidence of two active reconnection 
sites in the distant magnetotail, J. Geomag. 
Geoelectr., 48, p. 515, 1996. 

Huang, C.Y., C.K. Goertz, L.A. Frank, and G. Rostoker, 
Observational determination of the adiabatic index 
in the qiiiet time plasma sheet, Geophys. Res. Lett., 
16. pp. 563-566, 1989. 

Huang, C.Y., L.A. Frank, G. Rostoker, J. Fennell, and D.G. 
Mitchell, Nonadiabatic heating of the central plasma 
sheet at substorm onset, J. Geophys. Res., 97. pp. 
1481-1495,1992. 

Hughes, W.J., and D.G. Sibeck, On the three-dimensional 
structure of plasmoids, Geophys. Res. Lett., 14. pp. 
636-639, 1987. 

leda, A., S. Machida, T. Mukai, Y. Saito, T. Yamamoto, A. 
Nishida, T. Terasawa, and S. Kokubun, Statistical 
analysis on the plasmoid evolution with Geotail ob- 
servations, J. Geophys. Res., 103, p. 4,453, 1998. 

Jacquey, C., J.A. Sauvaud, and J. Dandouras, Tailward 
propagating cross-tail current disruption and dynam- 
ics of the near-Earth tail: A multi-point measurement 
analysis, Geophys. Res. Lett. 20(10), p. 983, 1993. 

Jorgensen, A.M., H.E. Spence, M.G. Henderson, G.D. 
Reeves, M. Sugiura, and T. Kamei, Global energetic 
neutral atom (ENA) measurements and their asso- 
ciation with the Dst index, Geophys. Res. Lett., 24, p. 
3173, 1997. 

Krimigis, S.M., and E.T. Sards, Energetic particle bursts 
in the Earth's magnetotail, in Dynamics of the Mag- 
netosphere, edited by S. Akasofu, pp. 599-630, D. 
Reidel, Dordrecht, Holland, 1979. 

Li, X., D.N.Baker, M.Temerin, G.D. Reeves, and R.D. 
Belian, Simulation of dispersionless injections and 
drift echoes of energetic electrons associated with 
substorms, Geophys. Res. Lett., 25(20), p. 3763, 1998. 

Lieberman, A., NASA-GSFC Magnetospheric Constel- 
lation Mission Document, NASA GSFC Code 730, 30 
March 1999a. 

Lieberman, A., Final Report: Nanosat Dispenser Ship 
Performance Feasibility Study for the Magneto- 
spheric Constellation Mission, NASA GSFC Code 
730, Swales Aerospace report SAI-RPT-284, 22 July, 
1999b. 

Lopez, R.E., and A.T.Y. Lui, A multi-satellite case study 
of the expansion of a substorm current wedge in the 
near-Earth magnetotail, J. Geophys, Res., 95, pp. 
8009-8017,1990. 

Lopez, R.E., H.E.J. Koskinen,T.I. Pulkkinen.T. Bosinger, 
R.W. McEntire, and T.A. Potemra, Simultaneous ob- 
servation of the poleward expansion of substorm elec- 
trojet activity and the tailward expansion of current 


sheet disruption in the near-Earth magnetotail, J. 
Geophys. Res,. 98(A6), p. 9285, 1993. 

Lui, A.T.Y., Plasma transport in the Earth's magnetotail, 
in Modeling Magnetospheric Plasma Processes, 
Geophys. Monogr. Ser., vol. 62, edited by G.R. Wilson, 
pp. 4 1 -54, AGU, Washington DC, 1991. 

Lui, A.T.Y, C.-L. Chang, A. Manofsky, H.-K. Wong, and D. 
Winske, A cross-field current instability for substorm 
expansions, J. Geophys. Res., 96, pp. 1 1389-11401, 1991. 

Lui, A.T.Y, Current disruption in the Earths magneto- 
sphere: Observations and models, J. Geophys. Res., 
101(A6), p.13067, 1996. 

Mauk, B.H., Quantitative modeling of the "Convection 
Surge” mechanism of ion acceleration, J. Geophys. Res., 
13, p.423, 1986. 

Maynard, N.C., W.J. Burke, E.M. Basinska, G.M. Erickson, 
W.J. Hughes, H.J. Singer, A.G. Yahnin, D.A. Hardy, and 
F.S. Mozer, Dynamics of the inner magnetosphere near 
times of substorm onsets, J. Geophys. Res., 101, pp. 
7705-7736, 1996. 

McComas, D.J., C.T. Russell, R.C. Elphic, and S.J. Bame, 
The near-Earth cross-tail current sheet: Detailed ISEE 
1 and 2 case studies, J. Geophys. Res., 91(A4). p. 4287, 
1 April 1986. 

Mcllvvain, C.E., Comments and speculations concerning the 
radiation belts, in Proceedings of the Joint IQSY/ 
COSPAR Symposium, London, 1 967, Part 1 , MIT Press, 
Cambridge, MA, USA, p.303, 1967. 

McPherron, R.L., C.T. Russell, and M.P. Aubrey, Satellite 
studies of magnetospheric substorms on August 15, 
1968, 9, Phenomenological model for substorms, J. 
Geophys. Res., 78(16), pp. 3131-3149, 1973. 

Miller, R.N., M. Ghil, and F. Gauthiez, Advanced data as 
simulation in strongly nonlinear dynamical systems, 
J. Atmos. Sci., 51, p. 1037, 1994. 

Moldwin, M.B., and W.J. Hughes, On the formation and 
evolution of plasmoids: A survey of ISEE 3 data, J. 
Geophys. Res., 97, p. 19259, 1992. 

Moore, T.E., R.L. Arnoldy, J. Feynman, and D.A. Hardy, 
Propagating substorm injection fronts, J. Geophys. Res. 
86 , p. 6713, 1981. 

Mukai, T., M. Fujimoto, M. Hoshino, S. Kokubun, S. 
Machida, K. Maezawa, A. Nishida, Y. Saito, T. Teresawa, 
and T. Yamamoto, Structure and kinetic properties of 
plasmoids and their boundary regions, J. Geomag. 
Geolectr., 48, p. 541, 1996. 

Nagai, T., Observed magnetic substorm signatures at syn- 
chronous altitudes, J. Geophys. Res., 87, pp. 4405-44 1 7, 
1982. 

Nagai, T., M. Fujimoto, Y. Saito, S. Machida, T. Teresawa, 

R. Nakamura, T. Yamamoto, T. Mukai, A. Nishida, and 

S. Kokubun, Structure and dynamics of magnetic 
reconnection for substorm onsets with Geotail obser- 
vations. J. Geophys. Res., 103, p. 4419, 1998. 

Nakamura, R., D.N. Baker, D.H. Fairfield, D.G. Mitchell, 
R.L. McPherron, and E.W. Hones Jr., Plasma flow and 
magnetic field characteristics near the midtail neutral 
sheet, J. Geophys. Res., 99, pp. 23591-23601, 1994. 


40 


Magnetospheric Constellation Mission... 


Nishida, A., T. Mukai, T. Yamamoto, Y. Saito, and S. 
Kokubun, GEOTAIL Observations of the reconnection 
process in the distant tail in geomagnetically active 
times, J. Geophys. Res., 100, p. 23663, 1995. 

Ohtani, S. p K. Takahashi. L.J. Zanetti, T.A. Potemra, R.W. 
McEntire, and T. lijima, Tail current disruption in the 
geosynchronous region, in Magnetospheric Substorms, 
Geophys. Monogr. Ser, vol. 64, edited by J. R. Kan, T.A. 
Potemra, S. Kokubun, and T. lijima, pp. 131-140, AGU, 
Washington DC, 1991. 

Ohtani, S., S. Kokubun, and C.T. Russell, Radial expan- 
sion of the tail current disruption during substorms: A 
new approach to the substorm onset region, J. Geophys. 
Res., 97, pp. 3129-3136, 1992. 

Ohtani, S. I., Earthward expansion of tail current disrup- 
tion: Dual -satellite study, J. Geophys. Res., 103 (A4). p. 
6815,1998. 

Owen, C.J., J.A. Slavin, I.G. Richardson, N. Murphy, and 
R.J. Hynds, Average motion, structure and orientation 
of the deep magnetotail determined from remote sens- 
ing of the edge of the plasma sheet boundary layer with 
E > 35 keV ions, J. Geophys. Res., 100, p. 185, 1995. 

Pulkkinen, T.I., et al., A study of magnetic field and cur- 
rent configurations in the magnetotail at time of a 
substorm onset, Plant. Space Sci., 39, p. 883, 1991. 

Quinn, J.M., and D.J. Southwood, Observations of 
parallelion energization in the equatorial region, J. 
Geophys. Res., 87, p. 10536, 1982. 

Reeves, G.D., R.W.H. Friedel, M.G. Henderson, A. Korth, 
PS. McLachlan, and R.D. Belian, Radial propagation 
of substorm injections, International Conference on 
Substorms-3, ESA SP-339, pp. 579-584, 1996. 

Robert, P, R. Gendrin, S. Perrault, and A. Roux, GEOS 2 
identification of rapidly moving current structures in 
the equatorial outer magnetosphere during substorms, 
J. Geophys. Res., 89, pp. 819-840, 1984. 

Roux, A., S. Perraut, A. Morane, P. Robert, A. Korth, G. 
Kremser, A. Pederson, R. Pellinen, and Z.Y. Pu, Role of 
the near Earth plasmasheet at substorms, In Magneto- 
spheric Substorms, Geophys. Monogr. Ser., vol. 64, ed- 
ited by J.R. Kan, T.A. Potemra, S. Kokubun, andT. lijima, 
pp. 201-214, AGU, Washington DC, 1991 . 

Schindler, K., A theory of the substorm mechanism, J. 
Geophys. Res., 79, p. 2803, 1974. 

Sergeev, V.A., T. Bosinger, and A.T.Y. Lui, Impulsive pro- 
cesses in the magnetotail during substorm expansion, 
J. Geophys. Res., 60, pp. 175-185, 1986a. 

Sergeev, V. A., R.J. Pellinen, T. Bosinger, W. Baumjohann, P. 
Stauning, and A.T.Y. Lui, Spatial and temporal charac- 
teristics of impulsive structure of magnetospheric 
substorm, J. Geophys. Res., 60, pp. 186-198, 1986b. 

Sergeev, V.A., R.J. Pellinen, and T.I. Pulkkinen, Steady 
magnetospheric convection: A review of recent results, 
J. Geophys. Res., 1996. 

Slavin, J.A., M.F. Smith, E.L. Mazur, D.N. Baker, T. lyemori, 
and E.W. Greenstadt, ISEE 3 observations of traveling 
compression regions in the Earths magnetotail , J. 
Geophys. Res., 98. p. 15425, 1993. 


Slavin, J.A., D.H. Fairfield, R.P Lepping. A. Szabo, M.J. 
Reiner, M. Kaiser, C.J. Owen, T. Phan, R. Lin, S. 
Kokubun, T. Mukai, T. Yamamoto, H.J. Singer.T. lyemori, 
and G. Rostoker, Simultaneous WIND, GEOTAIL and 
GOES 9 observations of magnetic field dipolarization 
and bursty bulk flows in the near-tail during a 
substorm, Geophys. Res. Lett., 24(8), p. 971, 1997. 

Slavin, J.A., D.H. Fairfield, M.M. Kuznetsova, C.J. Owen, 
R.P. Lepping, S. Taguchi, T. Mukai, Y. Saito, T. 
Yamamoto, S. Kokubun, A.T.Y. Lui, and G.D. Reeves, 
ISTP observations of plasmoid ejection: IMP 8 and 
Geotail, J. Geophys. Res., 103. p. 119, 1998. 

Song, Y., and R.L. Lysak, Turbulent generation of auroral 
currents and fields — A spectral simulation of two-di- 
mensional MHD turbulence, in Modeling Magneto- 
spheric Plasma, T.E. Moore and J.H. Waite (eds.),p. 197, 
AGU, Washington D.C., 1988. 

Spence, H.E., and M.G. Kivelson, Contributions of the low- 
latitude boundary layer to the finite width magneto- 
tail convection model, J. Geophys. Res., 98, pp. 15487- 
15496, 1993a. 

Spence, H.E., L. A. Frank, W.R. Paterson, and M.G. Kivelson, 
The Galileo Earth- 1 flyby: Comparison of the finite 
width tail convection model with observations, 1 993 Fall 
AGU meeting, EOS Trans. 74, p. 544, 1993b. 

Stern, D.P, and N.A. Tsyganenko, Uses and limitations of 
the Tsyganenko magnetic field models, 73, p. 489, 1992. 

Swift, D.W., Numerical simulation of the generation of elec- 
trostatic turbulence in the magnetotail, Geophys. Res., 
86, p. 2273, 1981. 

Taguchi, S., J.A. Slavin, M. Kiyohara, M. Nose, G.D. Reeves, 
and R.P. Lepping, Temporal relationships between 
midtail traveling compression regions and substorm 
onset: Evidence for near-Earth neutral line formation 
in the late growth phase, JGR, 103, 26607-26322, 1 998. 

Tajima.T., Astrophysical Plasmas, Frontiers in Physics se- 
ries, Addison Wesley, 1997. 

Talagrand, O., Assimilation of observations, an introduc- 
tion, J. Meteo. Soc. Japan, 75. p. 191, 1997. 

Tsyganenko, N.A., A magnetospheric magnetic field model 
with a warped tail current sheet. Planet. Space Sci., 
37. pp. 5-20, 1989. 

Tsyganenko, N.A., Modeling of twisted/warped magneto- 
spheric configurations using the general deformation 
method, J. Geophys. Res., 103 (A10), p. 23551, 1998. 

Tsyganenko, N.A., Modeling the Earth s magnetic field con- 
fined within a realistic magnetopause, J. Geophys. Res., 
100, p. 5599, 1995. 

Tsyganenko, N.A., and M. Peredo, Analytical models of the 
magnetic field of disk-shaped current sheets, J. Geophys. 
Res., 99, p. 199, 1994. 

Tsyganenko, N. A., and A.V. Usmanov, Determination of the 
magnetospheric current parameters and development 
of experimental geomagnetic field models based on data 
from IMP and HELIOS satellites, Planet. Space Sci., 
30, p. 985, 1982. 

Wang. X., A. Bhattacharjee. and A.T.Y. Lui, Collisionless 
tearing instability in magnetotail plasmas, J. Geophys. 
Res.,95(A9), p. 15047, 1 September 1990. 


Magnetospheric Constellation Mission... 


41 



Appendix B: Acronyms and Abbreviations 

MFE 

— 

Magnetic Field Experiment 




MHD 

— 

Magnetohydrodynamic 




MLI 

— 

Multi-Layer Insulation 

ACS 

— 

Attitude Control Subsystem 

mm 

— 

millimeter 

AFFCI/ 



MMS 

— 

Magnetospheric Multiscale 

CPS 

— 

A Formation Flying and Communications 

NASA 

— 

National Aeronautics and Space 



Instrument/Global Positioning System 



Administration 

AMIE 

— 

Assimilative Mapping of Ionospheric 

NENL 

— 

Near Earth Neutral Line 



Electrodynamics 

NRA 

— 

NASA Research Announcement 

ASIC 

— 

Application Specific Integrated Circuit 

NRO 

— 

National Research Opportunity 

BAI 

— 

Burst Accumulation Interval 

OP 

— 

Orbital Period 

BBF 

— 

Bursty Bulk Flows 

PSBL 

— 

Plasma-Sheet Boundary Layers 

BOL 

— 

Beginning of Life 

PVA 

— 

Plasma Velocity Analyzer 

CCE 

— 

Charge Composition Explorer 

QA 

— 

Quality Assurance 

CCSDS 

— 

Consulting Committee for Space Data 

RAI 

— 

Routine Accumulation Interval 



Systems 

RAM 

— 

Random Access Memory 

C&DH 

— 

Command and Data Handling 

RF 

— 

Radio Frequency 

CDPU 

— 

Command and Data Processing Unit 

RFP 

— 

Request For Proposal 

cm 

— 

centimeter 

rpm 

— 

Revolutions Per Minute 

CME 

— 

Coronal Mass Ejection 

sew 

— 

Substorm Current Wedge 

CPU 

— 

Central Processing Unit 

SEC 

— 

Sun-Earth Connection 

dB/K 

— 

decibel/Kelvin 

SECAS 

— 

Sun-Earth Connection/Advisory 

DMSP 

— 

Defense Meteorological Spacecraft 



Subcommittee 



Program 

SP 

— 

Spin Period 

DNL 

— 

Distant Neutral Line 

SST 

— 

Solid State Telescope 

DOD 

— 

Department of Defense 

STDT 

— 

Scientific and Technology Definition Team 

DPU 

— 

Data Processing Unit 

STI 

— 

Scientific and Technical Information 

DRACO 

— 

Dynamic Response and Coupling 

STP 

— 

Solar Terrestrial Probes 



Observatory 

TCS 

— 

Thermal Control System 

DS 

— 

Dispenser Ship 

UHF 

— 

Ultrahigh Frequency 

EIRP 

— 

Effective Isotropic Radiated Power 

UV 

— 

Ultraviolet 

EPA 

— 

Energetic Particles Analyzer 

W 

— 

Watt 

EPIC/ICS 

— 

Energetic Particles and Ion Composition/ 
Ion Composition System 




EPS 

— 

Electrical Power Subsystem 




ESA 

— 

Electrostatic Analyzer 




ESN 

— 

Essential Servicing Node 




FOV 

— 

Field of View 




FPGAs 

— 

Field Programmable Gate Arrays 




GaAs 

— 

Gallium Arsenide 




GEC 

— 

Geospace Electrodynamics Connection 




GFs 

— 

Geometric Factor 




g 

— 

gram 




GN2 

— 

Gaseous Nitrogen 




GN&C 

— 

Guidance, Navigation and Control 




GPS 

— 

Global Positioning Satellite 




GSFC 

— 

Goddard Space Flight Center 




HLBL 

— 

High-Latitude Boundary Layers 




HVPS 

— 

High-Voltage Power Supply 




IR 

— 

Infrared 




ISTP 

— 

International Solar Terrestrial Project 




Kb/sec 

— 

Kilobits/seconds 




LEP 

— 

Low Energy Plasma 




LLBL 

— 

Low-Latitude Boundary Layers 




m 

— 

meter 




MAG 

— 

Magnetometer 




MCM 

— 

Multi-Chip Module 




MCP 

— 

Micro-Channel Plate 




MCSTDT 

— 

Magnetospheric Constellation Science and 
Technology Definition Team 
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